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•'The  homicidal  chauffeur"  is  the  name  of  a  pursuit-evasion  differential  game 
originated  by  Isaacs  in  his  book.  Differential  Games.  In  this  game,  the  chauffeur  chases 
'a  slower  pedestrian  in  an  unbounded  parking  lot.  The  chauffeur's  control  is  his  turn 
rate,  bounded  in  magnitude,  and  the  pedestrian's  control  is  his  velocity  direction, 
which  can  be  changed  at  will.  The  pursuer  and  evader  seek  respectively  to  minimize  and 
maximize  the  capture  time,  when  the  radial  separation  becomes  less  than  a  known  capture 
radius.  The  two  equations  of  relative  motion  and  the  terminal  conditions  can  then  be 
written  in  terms  of  the  two  constant  parameters  of  the  game:  i)  the  speed  ratio,  ii)  tha 
ratio  of  capture  radius  to  pursuer's  minimum  turn  radius.  The  solution  to  the  problem 
(consists  in  finding  the  optimizing  strategies  of  both  players  as  functions  of  the 
position  relative  to  the  pursuer.  These  "min-max"  strategies  are  specified  in  terms 
of  the  local  position  variables  and  the  local  components  of  the  gradient  in  the  optimal; 
time-to-go,  values  of  which  are  known  at  the  termination  of  the  game.  The  equations  of* 
relative  motion  and  the  gradient  eqiations  are  then  integrated  retrogress ively  from 
termination,  so  as  to  fill  the  relative  space  with  optimal  paths  which  satisfy  all  of 
the  necessary  conditions  for  the  solution.  The  retrogressive  solution  is  made  dif- 
•ficult  by  the  presence  of  certain  "exceptional"  lines  in  the  relative  space.  These 
: lines  border  the  various  regular  regions  in  the  relative  space,  and  in  each  of  these 
‘regions  the  optimal  paths  take  a  specific  form.  The  exceptional  lines  themselves  may 
or  may  not  be  trajectories,  and  they  are  often  characterized  by  discontinuities  in  the 
{gradient  vector  and  in  one  or  both  of  the  controls.  In  addition  to  the  exceptional 
{lines  discovered  by  Isaacs  (e.g.,  the  "harrier",  across  which  the  time-to-go  is  dis¬ 
continuous,  the  "universal  line",  and  the  "equivocal  line") ,  two  new  types  of  lines  are 

«•"  jil-.n  "cuirrh  anw^l^p^"  and  "f^al  1  i  np " 
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L3.  Abstract  (continued) 

Furthermore ,  "safe-contact1*  motion,  for  which  the  evader  maneuvers  so  as  to  touch  but 
not  penetrate  the  capture  circle,  is  found  to  be  sometimes  optimal,  and  "dispersal 
lines"  for  both  pursuer  and  evader  are  encountered,  as  is  a  new  barrier.  Also,  cer¬ 
tain  details  in  the  preliminary  results  reported  at  the  1^  International  Differential  Game  (In¬ 
ference  in  1969  have  been  corrected.  Each  type  of  exceptional  line  (e.g.,  focal  line, 
evader's  dispersal  line,  etc.)  occurs  only  for  parameters  in  certain  regions  of  the 
parameter  space.  More  than  twenty  such  regions  are  found  to  exist,  each  of  which 
corresponds  to  a  different  configuration  of  exceptional  lines.  The  results  obtained 
are  believed  to  provide  a  complete  solution  to  the  game  for  all  combinations  of  the 
two  parameters.  \ 
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Chapter  I 


INTRODUCTION 


"The  homicidal  chauffeur"  is  a  pursuit-evasion  differential  game, 
first  described  by  Isaacs  in  his  epochal  book,  Differential  Games  (Ref. 

1).  As  one  of  the  first  physically  motivated  pioblems  in  this  emerging 
field  of  optimization,  it  has  both  historical  and  mathematical  interest. 
The  game  is  easily  described  and  visualized,  pitting  a  chauffeur  (the 
pursuer,  P  )  with  a  finite  minimum-turn  radius,  R  ,  against  a  slower 
pedestrian  (the  evader,  E  )  with  unrestricted  turn  radius.  The  pursuer's 
speed  advantage  is  partially  offset  by  the  evader's  agility;  these  fea¬ 
tures  provide  the  game  with  its  competitive  aspects. 

As  suggested  by  the  name  of  the  game,  P  wishes  to  minimize  the 
"capture  time,"  defined  as  the  time  when  the  mutual  separation  becomes 
less  than  the  capture  radius,  t  .  Simultaneously,  E  understandably 
seeks  to  maximize  this  time  or,  if  possible,  to  avoid  capture  forever. 

The  differential  game  thus  represents  a  generalization  of  the  minimum¬ 
time  optimization  problem  in  that  there  are  two  opposing  controls,  both 
of  which  Influence  the  relative  motion  and  the  resulting  optimal  or 
"min-max"  capture  time.  Much  of  the  interest  and  difficulty  of  the  game 
lies  in  the  fact  that  the  optimal  controls  are  not  always  unique,  and 
that  the  form  taken  by  the  optimal  solution  to  the  game  depends  very 
strongly  upon  the  speeds  of  the  players  and  the  maneuverability  of  P  . 
More  than  20  qualitatively  different  forms  of  the  solution  will  be  found 
to  be  possible,  these  forms  depending  on  the  relative  capabilities  of 
the  players.  The  unexpected  variety  and  complexity  of  the  solutions  to 
the  problem  support  and  Justify  certain  generalizations  regarding  dif¬ 
ferential  games  made  by  Isaacs  in  Ref.  1. 

1.1  Assumptions  and  Definitions 

The  game  is  visualized  as  taking  place  on  a  horizontal  plane.  When 
suitably  normalized,  the  differential  equations  giving  the  relative 
position  of  the  players  depend  on  two  parameters  which  are  the  speed 
ratio,  7  ,  and  the  ratio  of  capture  radius  to  P's  minimum  turn  radius, 


1 


(3  .  This  latter  parameter  is  used  in  defining  the  terminal  conditions  of 
the  game. 

An  Important  problem  in  a  pursuit-evasion  game  is  the  determination 
of  those  parameters  for  which  capture  is  possible  only  from  starting 
points  inside  a  certain  boundary,  the  evader  otherwise  avoiding  capture 
forever.  A  subsidiary  problem  arises  when  capture  is  known  to  be  possi¬ 
ble  from  all  starting  points,  but  for  which  there  exist  loci  of  discon¬ 
tinuity  in  the  optimal  time-to-go.  In  Chapters  II  and  III,  these  bar¬ 
riers"  are  shown  to  exist  for  certain  ranges  of  the  parameters. 

When  the  parameters  and  relative  position  are  such  that  capture  is 
possible,  a  solution  to  the  game  consists  in  specifying  the  optimal 
strategy  of  each  player  as  a  function  of  E's  position  (x,y)  relative 
to  P  .  These  "min-max"  strategies  are: 

i)  For  P  ,  the  rate  of  turn,  cp(x,y)  ,  with  |tp|  <  1  , 

li)  For  E  ,  the  velocity  direction,  V(x,y)  , 

and  when  capture  is  possible,  each  relative  position  is  associated  with 
an  optimal  positive  time-to-go,  V(x,y)  .  The  optimal  strategies  (which 
are  not  necessarily  unique)  are  such  that,  e.g.,  if  only  E  deviates 
from  his  optimal  strategy  for  a  finite  time,  capture  occurs  sooner  than 
V(x,y)  .  Likewise,  if  only  P  deviates  from  optimality,  capture  must 
occur  later  than  V(x,y)  .  This  is  the  source  of  the  "saddle-point" 
terminology  used  in  differential-game  theory. 

The  important  assumptions  of  the  game  are: 

i)  Both  P  and  E  know  the  relative  position  exactly. 

ii)  Both  players  know  the  speeds,  the  capture  radius,  and  P  s 
minimum  turn  radius. 

Other  assumptions  which  represent  departures  from  "reality"  (as  in,  say, 
air-to-air  combat  games)  are  evident  from  the  previous  discussion. 
Treating  the  speeds  as  constant,  forcing  the  motion  to  occur  in  a  plane, 
and  ignoring  any  response-time  lags  in  the  dynamics  of  P  and  E  are 
the  most  obvious  of  these. 
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It  is  appropriate  to  state  here  that  the  complexity  and  variety  of 
the  solutions  to  this  differential  game  were  underestimated  for  some 
time.  In  fact,  early  in  the  study  of  this  game,  it  was  mistakenly  be¬ 
lieved  that  little  remained  to  be  learned  from  it.  Attention  was  conse¬ 
quently  directed  to  generalizing  the  dynamical  model,  the  third-order 
"game  of  two  cars"  being  the  most  natural  generalization.  Certain  dif¬ 
ficulties  in  this  third-order  problem  raised  doubts  about  our  under¬ 
standing  of  the  simpler  homicidal  chauffeur  game.  These  doubts  proved 
to  be  well  founded,  and  led  to  the  decision  to  find  the  complete  solu¬ 
tion  to  the  present  game  for  all  values  of  the  parameters. 

The  pioneering  work  by  Isaacs  describes  a  treatment  of  the  problem 
which  is  complete  only  for  certain  rather  narrow  ranges  of  the  two  para¬ 
meters.  In  the  course  of  his  researches  into  this  problem,  various  loci 
in  the  relative  space  were  found,  across  which  either  or  both  of  the 
strategies  change.  Specifically,  a  barrier  and  an  equivocal  line  were 
found  to  enclose  "turn-away"  regions  in  the  relative  space,  and  the  y- 
axis  was  found  to  be  either  a  universal  line  or  a  pursuer's  dispersal 
line.  Loci  of  this  type,  which  separate  two  qualitatively  distinct 
families  of  relative  trajectories,  and  which  may  or  may  not  be  trajec¬ 
tories  themselves,  will  be  termed  "exceptional"  lines.  As  will  be  made 
clear  by  later  results,  Isaacs'  terminology  "singular  line"  has  an  un¬ 
desirable  connotation  to  most  control  theoreticians,  and  is  not  suffi¬ 
ciently  broad  in  meaning.  The  present  choice  is  made  to  accommodate  as 
large  a  class  of  such  arcs  as  possible,  and  its  definition  is  therefore 
intentionally  general. 

In  the  effort  to  extend  the  parameter  space  to  its  limits,  several 
new  types  of  exceptional  line  were  found  to  be  necessary.  Two  of  these 
lines  are  apparently  new,  not  only  to  the  specific  game  at  hand,  but  to 
the  theory  of  differential  games.  The  others  are  found  to  occur  in 
other  games  described  in  Ref.  1.  Each  type  of  exceptional  line  occurs 
only  for  parameters  in  certain  areas  of  the  parameter  space,  but  all 
such  areas  appear  to  have  been  found.  Consequently,  a  complete  solution 
to  the  game  is  believed  known,  for  all  values  of  the  parameters.  The 
startling  discoveries  made  at  intervals  during  this  study,  on  the  other 
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hand,  discourage  presumptions  of  this  type.  That  is,  there  may  yet  be 
novel  types  of  trajectories  or  strategies  which  might  occur  for  certain 
isolated  ranges  of  the  parameters.  So  many  cases  have  been  analyzed 
numerically  in  the  course  of  this  work,  however,  that  this  is  considered 
unlikely.  That  is,  the  deformations  in  the  trajectory  configurations 
are  found  to  be  smooth  and  continuous  with  changes  in  the  parameters,  so 
that  smooth  qualitative  changes  in  the  configurations  can  be  observed 
whenever  a  parameter  locus  is  crossed.  These  loci  are  actually  the  end 
result  of  the  study,  but  there  is  no  proof  that  the  loci  are  complete. 


1.2  Equations  of  Motion  and  General  Solutions 

The  pursuer's  location  in  the  playing  field  is  given  (Fig.  1.1)  by 
solving 


x  =  sin  $  =  b6  * 

p  P  P 

y  =  COS  a  =  C0  (1' 

Jp  P  P 

ep  ■  <p  .  M  £  1  • 

These  equations  effectively  normalize  the  dimensions  of  length  and  time 
as  P  has  been  given  unit  speed  and  a  maximum  turn  rate  of  unity.  The 
position  of  the  slower  evader  is  given  by  solving 


K  =  • 

where  0  <  7  <  1  •  1*  B's  control  is  taken  as  *  =  0.-6  p  .  the 

relative  coordinates  defined  as  in  Fig.  1.1  by 

x  =  (VXP)C*P  ‘  <ye“yP>Sep 

y  =  <vVseP  +  (Vyp)cep  * 


the  relative  position  <x,y)  is  given  by  the  solutions  to: 


*Here,  as  throughout  the  thesis, 
for  sine  and  cosine. 


s  and  c  are  used  as  abbreviations 
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x  =  -?y  +  ysy 

y  =  -l  +  <px  +  7cy  . 


(3) 


FIGURE  1.1.  P  and  E  in  Fixed  and  Relative 
Coordinate  Systems 


In  retrogressive  time,  dx  =  -dt  ,  and  the  equations  become 


*  =  cpy  -  7sy 

£  =  1  -  cpx  -  7cv  , 


(1) 


where  the  superscript  circle  denotes  the  derivative  with  respect  to  the 
dimensionless  time-to-go,  x  . 

When  capture  is  possible,  each  relative  position  (x,y)  has  asso¬ 
ciated  with  it  a  positive  constant,  the  optimal  time-to-go,  V(x,y)  , 
which  is  to  be  found  by  postulating  optimal  play  from  that  point  on. 

Its  total  time  rate  of  change  can  be  written  in  terms  of  its  gradients, 


5 


as 


dv  4  r„  ■ 

—  =  min  max  LV  x 

dt  cp  \|r  x 


+  V  y]  =  -1 

y 


(5) 


This  is  the  "main  equation"  of  Ref.  1,  and  using  (3)  it  may  be  expressed 
as  the  "Hamiltonian,"  or 

H  =  min  max  [-V  -  rp(V  y-V  x)  +  y(V  s\(r+V  c\|r)  ]  +  1=0.  (6) 

cp  y  x  y  x  y 


This  "min-max"  operation  gives  the  optimal  strategies  in  terms  of  the 

(known)  position  (x,y)  and  (unknown)  gradients  (V  ,V  )  , 

x  y 

(p  =  sgn  S  =  sgn(V  y-V  x) 
x  y 

(7) 

s±  _  cjL  _  1 

Vx  =  Vy  =  ' 


P  is  thus  to  turn  hard  left  or  hard  right,  unless  the  switch  func¬ 
tion  S  =  0  ,  and  E  is  to  run  normal  to  the  contours  of  constant  time- 
to-go.  Whenever  S  has  a  constant  sign  and  the  contours  of  constant  V 
have  continuous  second  partial  derivatives,  the  corresponding  region  is 
termed  "regular."  In  these  regions,  the  partial  derivatives  of  the  main 
equation  with  respect  to  x  and  y  must  vanish.  Thus,  for  example, 


a 


[V  x+V  y] 
x  y 


=  V  x  +  V  y 
xx  xy 


+  V  x  +  V  y 
xx  y  x 


=  0 


The  first  two  terms  on  the  right  equal  the  time  derivative  of  ,  and 
the  last  two  can  be  evaluated  by  reference  to  (3).  Repeating  this  opera¬ 
tion  with  respect  to  y  gives  the  gradient  equations  as 


V  =  -aH/dx  =  -mV 
x  y 


(8) 


vy  =  -aH/ay  =  tpvx  , 

where  cp  =  ±1  .  Note  that  since  H  is  stationary  with  respect  to  >)(  , 
the  partial  derivatives  are  most  easily  obtained  by  treating  t  >  as  well 

as  CD  and  (V  ,V  )  ,  as  fixed. 

Y  x  y 
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In  retrogressive  time,  and  for  general  initial  conditions,  the  solu¬ 
tions  to  (8)  are 


VX  =  j  'l  +  Vy  S  C+ot<pT> 

Vo 

Vy  =  Jt  +  Vy  C<V*T>  • 
Vo  o 


(9) 


We  have  found  that,  when  S  >  0  or  S  <  0  ,  E's  optimal  direction 

in  relative  space  varies  linearly  with  time,  -f  =  +  pT  .  Since  P's 

o 

rate  of  turn  is  ,  this  means  that  E's  motion  is  rectilinear  in  real 
space,  and  that  the  linear  time  variation  is  due  to  P's  turn  strategy. 

The  retrogressive  solutions  to  the  state  equations  (4),  over  any 
time  interval  during  which  p  is  constant  (  S  of  constant  sign)  are 


x  =  cp(l-cx)  +  x^ct  +  py^st  -  yts (ij/^+cpt) 
y  =  (l-q>xo)sT  +  y^ct  -  ytc (ij/^+cpt)  , 


(10) 


where  the  controls  \|r  and  tp  are  given  by  (7)  in  terms  of  the  gra¬ 
dients  (vx0>vy0)  •  Determination  of  the  gradients,  as  functions  of  the 
relative  position  (x,y)  ,  is  the  most  important  source  of  theoretical 
difficulties  in  the  problem,  since  these  functions  used  in  (7)  essen¬ 
tially  yield  the  feedback-solution  to  the  problem. 


The  approach  taken  in  the  investigation  of  the  game  is  the  same  as 
that  used  by  Isaacs  in  Ref.  1.  That  is,  terminal  conditions  for  x  and 
W  are  found,  in  terms  of  the  parameters  '  and  7  ,  which  lead  to 
terminal  values  for  the  controls  \|i  and  p  .  Then  the  retrograde  solu¬ 
tions  (10)  are  used  to  find  where  these  terminal  points  must  have  come 
from.  Since  any  point  outside  the  capture  circle  is  a  possible  initial 
condition,  the  intention  is  simply  that  of  filling  the  entire  space  with 
retrogressive  paths.  We  will  lind,  however,  that  for  any  choice  of  para¬ 
meters,  the  space  is  subdivided  by  exceptional  lines  into  various  regular 
regions,  in  each  of  which  (10)  is  satisfied.  The  exceptional  lines 
which  border  the  regular  regions  have  one  or  more  of  the  following  dis¬ 
tinguishing  features: 
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i)  The  time-to-go  is  discontinuous  across  the  line; 

ii)  The  gradient  in  the  time-to-go  is  discontinuous  across  the 
line  (as  is  therefore  E's  strategy); 

iii)  The  switch  function  on  the  line  has  the  value  S  =  0  ; 

iv)  The  line  is  the  capture  circle  itself. 

These  lines  may  or  may  not  be  optimal  paths  themselves.  Certain  addi¬ 
tional  necessary  conditions  must  be  satisfied  by  the  exceptional  lines 
which,  together  with  the  main  equation,  will  guarantee  global  optimality. 
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Chapter  II 


CLASSICAL  RESULTS 

Isaacs'  analysis  of  this  game,  as  given  in  Ref.  1,  begins  with  the 
determination  of  the  "usable  part"  of  the  capture  circle.  This  allows 
the  specification  of  the  strategies  (<p,^)  and  the  terminal  conditions 
(xo.yo)  in  the  solutions  given  by  (10)  of  the  previous  section.  Work¬ 
ing  retrogressively  from  termination,  he  then  shows  how  to  fill  the 
(x,y)  space  with  trajectories,  using  the  main  equation  and  its  deriva¬ 
tives.  Implicit  in  this  analysis  is  the  restriction  of  parameters  to 
certain  regions  in  the  parameter  space,  which  will  be  specified  below. 
But  these  "classical"  results,  apart  from  their  many  interesting  fea¬ 
tures,  provide  a  foundation  of  conclusions  which  will  ultimately  lead  to 
solutions  in  adjacent  areas  of  parameter  space. 

2 . 1  One-Stage  Game 

The  "usable  part"  of  the  capture  circle,  on  which  the  optimally 
played  game  must  end,  is  easily  determined  by  writing  the  main  equation 
in  polar  coordinates, 

min  max[V  r  +  V  0]  =  -l  m 

q>  r  p 


where 

r  =  — c  $  +  jc  (y-e) 

6  =  -Ur  +  [s0  *  7s(\|/-0)]/r 

Since  termination  is  defined  independent  of  the  terminal  angle  $ 

o 

the  terminal  value  of  the  angular  adjoint  V  must  be  zero,  while  the 

0 

radial  adjoint  at  this  time  is  given  by  substituting  for  r  in  (1)  and 
performing  the  indicated  maximization: 

max  V  [-c8  +  7cty-0  )]  =  V  (7-c0  )  =  -1  .  (2) 

^  r  o  T  o  r  o 

This  brief  analysis  has  provided  three  important  results: 
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i)  At  termination,  when  r  =  p  ,  P's  control  (turn  rate)  is 
irrelevant,  and  E's  control  (direction)  is  radially  out¬ 
ward,  \|/  -  ,  as  suggested  by  intuition; 

ii)  The  terminal  udjoints  are  V  =  0  ,  V  =  — .  >  0  ; 

6  t  c0o-/ 

iii)  The  terminal  value  of  the  angular  coordinate  must  satisfy 

Iri  I  <  8  ,  where  the  "usable  part"  is  defined  by  the  half- 

UP  A  -1 

angle  0  =  cos  y  . 

up 


At  the  terminal  location  r  =  £  ,  9  =  9q  ,  we  first  find  the  switch 
function  which  determines  P's  strategy,  according  to  (7)  of  Sec.  1.2. 
The  cartesian  adjoints  are 


V 

x 


V  80 

r 


g  Q 

+  <V0/r)cS  "  ^ 


V  =  V  c0 
y  r 


(V  /r?s0  = 
& 


(3) 


With  these  terminal  values  for  the  gradients,  it  is  easy  to  verify  that 
the  switch  function  is  the  angular  gradient,  and  therefore  it  vanishes. 
The  retrogressive  time-derivative,  however,  is 


80 

V  ’ 


(4) 


which  means  that  y  -  sgn  S  =  1  ,  for  >  0  ;  hence  P  is  turning 
toward  E  just  before  termination.  This  result  has  an  obvious  intui¬ 
tive  interpretation.  More  generally,  we  note  that,  according  to  (7)  of 
Sec.  1.2,  whenever  S  =  0  ,  E's  velocity  is  radial;  E  is  running 
directly  toward  or  away  from  P  . 

For  right  half-plane  terminal  conditions,  xq  =  (5s0q  >  0  ,  the 
solution  (10)  of  Sec.  1.2  is,  with  0=1, 

x  =  1  -  ct  +  O-7t)s(0  +t) 

°  (5 

y  =  st  +  (S-7t)c(0o+t)  . 
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Here  it  is  noted  that,  when  t  =  3/7  ,  the  point  (x,y)  is  inde¬ 
pendent  of  the  arrival  angle  9q  ,  and  that  this  point,  A  ,  has  coor¬ 
dinates 


x 


A 


1 


-  CT. 


y 


A 


ST. 


(6) 


where 


ta  -  • 

These  equations  imply  that  the  point  A  is  located  on  a  unit  circle 
in  the  x,y  plane  according  to  the  value  taken  by  the  ratio  3/7  ,  as 
shown  in  Fig.  2.1. 


FIGURE  2.1.  Locus  of  Point  A  for  3>  =  +1 
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We  note  here  that  all  optimal  trajectories  terminating  at  6  <6 

o  —  ue 

tiln  Lie  considered  as  emanating  from  this  point  A  ,  which  is  then  a 
dispersal  point  ior  E  .  This  is,  in  fact,  the  convention  adopted  by 
Isaacs,  and  according  to  the  arrival  angle  "preferred"  by  E  ,  his  op¬ 
timal  strategy  \|/^  at  the  dispersal  point  need  not  be  unique.  This 
strategy  need  only  satisfy  the  inequality, 

3/7  <  VA  <  3/7  +  cos-1 7  ,  (7) 

any  one  of  which  leads  to  r  =  3  in  the  time  3/7  . 

As  shown  in  (1)  when  r  =  3  ,  the  radial  velocity  is  r  =  7  -  c0 

for  the  paths  V-' .  and  hence  r  =  0  when  8=6  .  The  sidemost  path 

o  up 

from  A  ,  as  shown  in  Fig.  2.2(a),  touches  but  does  not  penetrate  the 
capture  crclt  Accordingly,  this  path  is  here  considered  as  distinct 
from  its  11..  t  ghnof' ,  in  that  it  leads  to  "safe-contact"  for  E  ,  rather 
than  capture.  Thus,  our  convention  differs  from  that  of  Isaacs  in  that 
we  consider  E  to  have  only  one  strategy  from  A  ,  assuming  that  P 
turns  toward  E  ,  and  this  will  be  called  the  "barrier"  or  "safe-contact" 
strategy. 

It  is  also  interesting  to  see  that,  if  P  turns  toward  E  when  he 
is  initially  on  IB  but  above  the  point  A  ,  E  chooses  the  barrier 
strategy.  This  path  arrives  earlier  at  the  capture  circle  than  does  the 
two-stage  trajectory  via  the  y-axis,  but  here  it  is  simply  evident  that 
E  prefers  an  early  safe-contact  to  a  later  arrival  with  penetration, 
and  he  runs  accordingly.  Whether  or  not  P  also  prefers  the  barrier 
strategy  will  be  determined  by  a  certain  parameter  inequality,  to  be 
derived  in  Sec.  2.3. 
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f 


2.2  Universal  Line 


In  Fig.  2.2,  trajectories  from  the  neighborhood  of  A  terminate  on 
the  capture  circle  over  a  range  of  angles  for  cp  =  1  .  A  symmetric  group 
of  paths  for  which  cp  =  -1  arrives  at  the  capture  circle  from  the  left. 
One's  curiosity  is  therefore  aroused  regarding  the  conditions  under 
which  cp  =  0  is  optimal.  We  have  shown  that  cp  =  sgn  S  ,  and  the  clas¬ 
sical  definition  of  universal  line  (UL)  or  "singular  arc,"  is  a  path 
along  which  the  switch  function  is  identically  zero.  Thus,  using  the 
equations  of  motion  (3)  and  adjoint  equations  (8)  of  Sec.  1.2,  we  find 


S=Vy-Vx  =  0 
x  y 


(1) 


s  =  v^d-ycijc)  +  v^ysy  =  0  . 

In  order  for  V  and  V  not  to  be  both  zero,  these  equations  imply 
x  y 

that  the  determinant  of  their  coefficients  is  zero,  or 

x  -  y  (xci|r-ys\|/)  =  0  .  (2) 

But,  since  we  have  shown  in  (7)  of  Sec.  1.2  that  E's  control  (s<|r,c\|i) 

is  proportional  to  the  gradient  W  =  (V  ,V  )  ,  the  parenthesis  in  (2) 

x  y 

is  proportional  to  S  ,  and  therefore  it  vanishes,  leaving 

x  =  0  (3) 


as  the  only  candidate  for  a  UL  in  the  field  of  primaries.  Thus,  the 

hypotheses  in  (1)  are  fulfilled  only  for  motion  of  E  down  the  y-axis. 

On  this  line,  the  switch  function  (1)  Implies  cp  =  V^y  =  0  ,  and  hence 

V  =  tc  =  0  ,  so  that  both  P  and  E  follow  the  same  straight  path  in 
x 

real  space. 

The  condition  (3)  is  necessary  but  not  sufficient  for  the  existence 
of  the  UL  ,  and  in  fact  w6  will  find  that  portions  of  the  y-axis  can 
instead  be  a  dispersal  line  for  F  ,  or  a  dispersal  line  for  E  ,  de¬ 
pending  again  on  6  and  7  •  For  parameters  such  that  a  barrier  exists 
(to  be  specified  in  Sec.  2.3),  the  trajectory  configuration  near  the 
positive  y-axis  resembles  the  upper  paths  of  Fig.  2.2.  As  shown  in 
Ref.  1  (p.  194),  these  tributary  paths  correspond  to  sharp  turns  by  P  , 
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with  E  running  tangent  to  P's  turn  circle.  The  contours  of  time-to- 
go,  shown  as  dashed  lines  in  the  figure,  are  smooth  at  the  y-axis,  which 
implies  that  E's  strategy  is  continuous  here,  being  the  same  straight- 
line  motion  before  and  after  P  stops  turning. 

Specifically,  for  a  trajectory  which  encounters  the  y-axis  at  the 
ordinate  y^  ,  the  retrogressive  continuous  solutions  to  the  adjoint 
equations  are,  with  cp  =  +1  , 


V 

x 


V 

y 


G) 


and  the  path  equations  are 

x  =  1  -  ct  +  (y1-7t)st 


(5) 


y  =  st  +  (y1-7t)cr 

where  T  is  measured  on  the  path  back  from  the  ordinate  yx  .  The 
magnitude  of  the  gradient  vector  along  and  adjacent  to  the  UL  is  given 
by  the  main  equation  when  x  =  =  0  ,  so  |w|  =  Vy  =  =  I/O”?)  • 

For  the  rectilinear  chase  to  which  this  result  applies,  the  interpreta¬ 
tion  of  the  gradient  is  particularly  simple;  i.e.,  the  change  in  optimal 
time-to-go  resulting  from  a  change  in  relative  (y-axis)  position  is  the 
inverse  of  the  difference  in  speeds. 


2.3  Barrier 

When  the  trajectories  derived  in  (5)  of  Sec.  2.1  have  the  boundary 

value  0=0  .the  resulting  locus  defines  the  "barrier,"  S  .  There 
up 

is  here  an  infinite  discontinuity  in  the  radial  adjoint,  as  shown  by- 
putting  0=6  =  cos"1?  in  (2)  of  Sec.  2.1.  At  the  barrier,  then, 

according  to  (3)Pof  Sec.  2.1,  both  Vx  and  Vy  are  infinite  in  magni¬ 
tude.  Accordingly,  as  the  barrier  is  approached,  the  main  equation  may¬ 
be  divided  by  the  magnitude  of  the  adjoint,  so  that,  with  Vg  ~  0  ’ 
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lim 

V  -»  co 
r 


min  max 

<p  v 


[V  x+V  y] 
x  yJJ 


lim 

V  -♦  00 
r 


HL_ 

V 

r 


=  0  . 


Hence,  defining  a  normalized  adjoint  vector, 


Vv  =  (v  ,  V  )  = 


v/vf 


(V  ,v  )  . 


,+Vf, 


(1) 


the  main  equation  becomes 


min  max[v  x+v  y]  =  0 

cp  \jf  x  y 


(2) 


The  adjoint  vector  'v'v  is  therefore  normal  to  the  barrier  trajectory, 
as  is  E's  optimal  velocity  direction,  \|i  . 

Thus,  for  r  =  p  ,  0Q  =  6up  ,  (2)  becomes,  with  (vx,Vy)  = 
(s0up’CSup)  and  ^  =  +1  • 


max[s0  (Pee  -7s*  )  +  c0  (1-BsS  -?cti  )]  =  0  (3) 

^  up  ^  up  To  up  K  up  Yo  J 

o 

which  implies  =  G 

o  up 

The  retrogressive  adjoint  equations  are,  with  their  boundary  condi¬ 
tions, 


V  =  V 
x  y 


V  =  -v 
y  x 


v  (o)  =  s e 
x  up 


Vy(0)  =  00  up  • 


(4) 


and  the  solutions  give  E's  barrier  strategy  for  cp  =  +1  , 

v  =  s(g  +t) 
x  up 


Vy  "  C(eup+T)  - 


so  that 


(5) 


*( T)  =  ®up  +  T 


(6) 


The  equations  of  the  barrier  are  then  given  by  (5)  of  Sec.  2.1, 
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(7) 


x  =  1  -  ct  +  (p-7i)s(eup+T) 

y  =  st  +  (p-7t)c  (6up+T)  . 

The  barrier  is  termed  a  semipermeable  surface  by  Isaacs  (p.  70)  since, 
if  cp  =  +1  i  no  strategy  of  E  allows  x  to  penetrate  the  barrier  in 
the  "E-direction,"  and  if  \Jr  is  given  by  (6),  no  strategy  of  P  allows 
penetration  in  the  P-direction.  Here  the  P-direction  refers  to  the  side 
of  smaller  time-to-go,  and  in  general  the  optimal  time-to-go  is  discon¬ 
tinuous  at  the  barrier.  Thus,  for  x  on  the  barrier,  a  departure  from 
the  strategy  \|r  =  0up  +  T  by  E  results  in  motion  away  from  the  bar¬ 
rier  in  the  P-direction,  to  a  region  of  small  V  .  Likewise,  a  depar¬ 
ture  by  P  from  cp  =  +1  ,  if  E  is  paying  attention,  moves  E  across 
the  barrier  to  a  region  of  large  (perhaps  infinite)  V  ,  The  discon¬ 
tinuity  in  V  decreases  to  zero  at  the  end  B  of  $  ,  which  is  indi¬ 
cated  by  the  vanishing  of  P's  switch  function.  Combining  (5)  and  (7), 
this  function  is 

S  =  vy-vx  =  7-  c(6  +t)  .  (8> 

x  y  up 

P's  barrier  strategy  is  thus  rp  =  +1  for  S  >  0  ,  or  for 

0  <  T  <  2(n-eup)  £  tB  ,  (9) 

and  unless  the  barriers  intersect  on  the  positive  y-axis,  T0  is  the 
time  required  to  traverse  the  barrier. 

The  right  and  left  barriers  Intersect  tangentially,  on  the  y-axis, 
for  given  values  of  p  and  7  ,  if  T  exists  such  that  x(j3,7if)  — 

8(P,  7, t)  =  0  .  Using  (5)  of  Sec.  2.1,  and  eliminating  T  from  these 
equations,  gives  the  condition  for  tangential  closing  as 

P  =  7  sin  *7  +  (10) 

Implicit  in  this  derivation  is  an  important  parameter  inequality, 
necessary  for  the  existence  of  the  barrier.  This  inequality  can  be 
derived  by  calculating  the  velocity  components  when  Cp  =  +1  and 
q  =  COB'17  at  r  =  P  .  In  polar  coordinates,  these  velocities  are,  in 
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general , 


r  =  -cO  +  7c  (ty-0) 

e  =  -<p  +  i[se  +  7s(^-e)  ] 


(XI) 


and  for  the  given  terminal  conditions,  the  retrogressive  velocity  com¬ 
ponents  are  found  to  be 


r  =  0 


O 

e 


-l 


(12) 


The  retrogressive  radial  acceleration,  however,  is 


V  =  - 9 

and  will  be  positive  only  so  long  as  0  <  0  or,  by  (12),  when 

p2  +  y2  <  1  •  (13) 

This  is  the  parameter  space  inequality  necessary  for  the  existence  of 
the  barrier  which  terminates  at  0  =  0up  with  tangential  velocity. 

For  p  larger  than  the  value  given  in  (10),  the  barriers  do  not 

contact  the  y-axis  but  instead  curve  away  to  the  rear,  according  to  the 

parametric  equations  (7).  In  this  case,  the  end  of  the  barrier  is 

reached  at  =  2(n-cos-1?)  ,  and  this  point  has  coordinates  given  by 
B 

(7);  i.e., 


*B  [2  ^  -  P  +  7Tb] 

yB  =  -7[2  -  e  +  7Tb]  • 


(14) 


Note  that  the  bracket  in  (14)  is  the  radial  distance,  rg  ,  and  that 
J(r  -9)  =  Vl-72  +  7(rt-cos  1y)  -  3  >  0  . 

D 


This  inequality  is  implied  by  (13),  and  it  means  that  the  end  of  the 
barrier,  B  ,  is  always  outside  the  capture  circle. 
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It  is  possible  for  parameters  to  be  such  that  the  point  A  is  not 
on  the  barrier;  i.e.,  >  Tg  .  Substituting  for  these  retrogressive 

times  from  (9)  and  from  (6)  of  Sec.  2.1,  this  inequality  can  be  reduced 
to 

c  (3/27)  +  7  >  0  .  (15) 


2*4  Equivocal  Line  and  Pursuer's  Dispersal  Line 

For  parameters  such  that  the  barriers  are  open,  the  equivocal  line 
(EL)  joins  the  end  of  the  barrier  to  the  y-axis  or  to  the  capture  circle, 
and  the  y-axis  below  the  EL  is  a  dispersal  line  for  P  .  The  EL  is 
a  locus  of  points  along  which  E's  strategy  is  locally  the  same  as  his 
incoming  strategy,  which  is  such  as  to  cause  P's  switch  function  to  be 
zero;  i.e., 


V  y  -  V  x  =  0  . 
*  y 


(1) 


This  implies  that  tan  ip  =  x/y  ,  and  that  E  is  running  directly 
toward  P  ,  Though  the  minimization  through  the  main  equation  yields 
cp  =  sgn  S  ,  the  signum  operator  is  indeterminate  at  S  =  0  ,  and  so  we 
require  the  intermediate  value  of  turn  rate  by  other  means.  Continuity 
of  the  components  of  the  gradients  along  the  EL  ,  and  requiring  unit 
rate  of  decrease  of  the  time— to-go  thereon  together  gives 

VeL  +  V/EL  =  V^EL  +  Vy*EL  =  '  (2) 


These  equations  are  interpreted  vectorially  in  Fig.  2.3,  where  it  is  seen 
that  the  adjoint  vectors  before  and  after  arrival  at  the  EL  hove  equal 
components  along  the  EL  ,  and  that  the  main  equation  is  satisfied  both 
above  and  below  the  EL  . 


Thus,  the  EL  is  a  locus  along  which  E  can  choose  between  two 
strategies,  and  P's  strategy  depends  on  this  choice.  If  E  chooses 
^  ,  P  is  forced  to  mix*  his  strategies,  for  a  resultant  intermediate 
control,  cp(x,y)  ,  depending  on  E's  position.  This  "chatter"  strategy 


*The  "mix"  is  accomplished  at  infinite  frequency,  such  that  p  =  2p-l  , 
where  p  =  proportion  of  time  spent  at  ?  =  +1  . 
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FIGURE  2.3.  Velocities  and  Gradients  Near  the  EL 


produces  a  path  for  P  of  time-varying  curvature  in  real  space.  When 
E  chooses  \|f+  ,  he  flees  along  the  tangent  to  P's  minimum  turn  circle, 
and  P's  optimal  strategy  jumps  to  <p  =  +1  .  For  EL  motion,  the 
retrogressive  equations  read 


%  =  yy  -  7sy 

$  =  1  -  $x  -  7c\|r 


(3) 


where 


(4) 


The  main  equation  on  the  prior  side  then  gives  the  magnitude  of  the 
adjoint  vector,  so  that,  with  (4), 
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(5) 


V 

x 


x 

y+7r 


V  =  — . 

y  y+7r 

To  evaluate  Cj5  ,  (4)  is  used  to  give  E's  control  in  (2),  which 
then  becomes 


v^(<py+yx/r)  +  V  (l-$x+7y/r)  =  1  , 


(0) 


where  the  gradients  are  to  be  found  as  functions  of  (x,v)  on  the  EL  . 
Specifically,  the  adjoints  and  position  are  given  by 


V 


+ 

x 


ST 


1 


1-7 


x  =  1  -  ct  +  (y1~7T1)sx1 
y  =  stj  +  (y1-7t1)cT1  , 


(7) 


(8) 


where  T^  is  measured  back  from  the  position  (0,y^)  .  Eliminating 
(y  1» 1 )  from  (8)  and  solving  the  resulting  quadratic  equation  give 
the  time  from  EL  to  y-axis  as  a  function  of  position  on  the  EL  : 


1-x+hy 
c  T  =  1 

1  v.2  ! 

h  +1 


y-h(l-x) 

h2+l 


where  h2  =  x2+y2-2x  .  Substituting  these  relations  into  (7)  and  this 
result  into  (6)  give  P's  intermediate  control  as 


$  =  1 


2 

r  -x+hr 

7  TTr^yT  * 


(9) 


It  is  easy  to  prove  that  $  t  1  ,  since  the  numerator  in  (9)  cannot 
vanish  except  at  x  =  y  =  0  .  At  the  left  end  of  the  EL  ,  x  =  0  ,  and 
the  chatter  strategy  here  is  C?  =  1-7  •  Typical  numerical  examples  pro- 
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duue  deviations  of  cp  over  the  EL  of  only  a  f  ew  percent  from  this 
initial  value.  The  point  B  ,  where  the  EL  joins  the  barrier,  has 
coordinates  given  by  (14)  of  Sec.  2.3.  These  are  the  initial  conditions 
for  the  retiogressive  integration  of  the  EL  ,  which  then  proceeds  by 
using  (4)  and  (9)  in  the  equations  (3). 

It  can  readily  be  shown  that  the  slopes  of  EL  and  barrier  are 
equal  at  B  .  The  reason  is  that  y  +  7r  =  0  at  the  point  B  ,  because 

r  =  0  on  both  arcs  and  therefore  the  slope  at  B  is  independent  of  cp  , 

dy/dx  = v 1 -T2/?  •  This  also  implies,  according  to  (5),  that  the  gradient 
becomes  unbounded  at  B  ,  since  the  sensitivity  of  capture  time  to  a 
change  in  position  is  infinite  at  this  point# 

Our  convention  regarding  the  barrier  has  been  briefly  discussed  in 
Sec.  2.1.  While  this  topic  is  discussed  completely  in  Chapter  III, 
special  consideration  must  be  given  here  to  its  end-point,  B  .  For 
those  trajectories  which  begin  above  the  EL  ,  we  have  found  that  E  can 

choose  to  remain  on  the  EL  to  the  point  B  ,  if  P  plays  optimally. 

But,  since  E  can  avoid  termination  by  playing  the  barrier  strategy,  P 
must  prevent  E  from  reaching  B  .  As  mentioned  by  Isaacs  (p.  300),  P 
can  postpone  his  switch  to  q5  by  an  infinitesimal  time.  The  ensuing 
equivocal  path  is  then  imagined  as  occurring  slightly  below  the  EL  , 
such  that  B  cannot  be  reached  by  E  . 

As  shown  in  Fig.  2.5(b),  a  portion  of  the  negative  y-axis  may  be  a 
UL  for  P  ,  With  cp  =  0  .  That  is,  if  the  EL  encounters  the  negative 
y-axis  (instead  of  the  capture  circle),  the  trajectories  leading  to  the 
EL  from  the  "turn-away"  region  arrive  also  at  the  negative  y-axis  just 
under  the  capture  circle.  This  segment  will  then  be  a  UL  as  discussed 
in  Bee.  2.2,  and  this  "reverse  chase"  strategy  of  E's  is  continuous 
with  the  strategies  on  either  side  of  the  y-axis,  as  E  strives  to 
delay  P’s  turn  maneuver. 

The  intersection  of  the  UL  and  EL  is  a  multiple-choice  point  in 
the  relative  space.  That  is,  both  P  and  E  must  choose  strategies, 
and  trajectories  can  depart  from  this  point  in  four  different  directions. 
If  E  chooses  to  continue  his  reverse  chase  policy  at  the  EL  ,  the 


22 


switch  function  remains  at  zero  while  P  executes  a  slow  turn  at  the 
intermediate  rate,  $(x,y)  .  If  he  does  not  so  choose,  E  runs  tangent 
to  the  minimum  turn  circle  chosen  by  P  •  On  arrival  at  this  four-way 
junction,  then,  P  "first"  chooses  a  turn  direction  (or  the  sign  of  <p  ) , 
and  K  ,  by  his  choice  of  strategies,  simultaneously  forces  P  to  turn 
at  one  of  two  specific  magnitudes  (|cp|  =  <p  or  1)  . 

Immediately  below  the  intersection  of  the  EL  with  the  y-axis 
(where  y  is  given  by  numerical  solution  of  the  EL  equations),  we 
find  the  value  of  the  switch  function  to  be  S+  =  V^y  ,  Using  (7)  and 
(8),  this  reduces  to  S+  =  (l-CT^/Q-y)  .  Substituting  for  c^  in 
terms  of  the  position  (0,y)  then  yields 

s+  =  _ -  >  o  .  (1°) 

(1-7)  d+y2) 

A  similar  analysis  for  trajectories  in  the  left  half  plane  (<p  =  -1) 
shows  that  S~  <  0  ,  and  since  this  analysis  can  hold  for  all  retrograde 
paths  not  intersecting  the  EL  ,  we  conclude  that  the  switch  function  is 
discontinuous  across  the  negative  y-axis  at  all  points  below  the  EL  . 

A  dispersal  line  is  defined  in  Ref.  1  as  a  locus  of  initial  condi¬ 
tions  along  which  either  or  both  of  the  strategies  are  not  unique.  They 
are  often  found  as  the  retrograde  intersection  of  two  distinct  families 
of  optimal  paths,  for  each  of  which  the  main  equation  is  satisfied.  At 
the  intersection,  the  optimal  time-to-go  is  the  same  for  either  pair  of 
strategies,  and  the  dispersal  line  condition  is  expressed  in  terms  of 
the  two  fields  of  trajectories,  as  follows: 

W“.x*  =  W+.x+  .  (11) 

This  reduces  to  a  differential  equation  for  the  slope  of  the  dispersal 
line, 

v"-v+ 

dy  _x x  (12 ) 

~T~  =  ~  ,  > 

d*  v--V 

y  y 

where  the  superscripts  refer  to  opposite  sides  of  the  line. 
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The  discontinuity  of  S  across  this  portion  of  the  negative  y-axis, 
as  derived  in  (10),  coupled  with  the  symmetry  of  strategies  and  paths 
across  it,  is  a  special  case  of  (12),  as  here  we  have  Vx  =  -V^  and 
V-  =  v+  This  portion  of  the  y-axis  is  then  a  dispersal  line  along 

which  the  sign  of  S  is  undefined.*  A  diagram  of  the  velocity  vector 
components  for  x  on  the  negative  y-axis  below  the  EL  shows  that 
P's  contribution  to  the  x-velocity  (— Cpy )  exceeds  E's  contribution 
(ys\|f)  ,  and  hence  P  determines  the  departure  direction.  The  line  is 
therefore  called  a  "pursuer's  dispersal  line"  (PDL)  .  The  present  case 
has  a  geometrically  obvious  interpretation  (when  E  is  directly  behind 
P  ,  P  can  turn  sharp  right  or  sharp  left,  and  E’s  strategy  depends  on 
P's  choice)  but  we  will  find  other  more  general  circumstances  in  which 
(12)  holds.  In  these  other  situations,  the  discontinuity  W  ^  W 
implies  a  discontinuity  in  E's  control  4+  ^  V  >  »  but  11  need  not 
imply  a  change  in  P's  control.  That  is,  the  switch  function  can  be 
discontinuous  across  the  line,  while  retaining  the  same  sign  on  either 
side. 


2.5  Parameter  Space  and  Trajectories 

To  summarize  the  developments  of  this  chapter,  it  has  been  shown 
that  the  "usable  part"  of  the  capture  circle  is  an  angular  sector  on 
which  the  game  must  end.  If  the  parameters  are  such  that  a  trajectory 
arrives  tangentially  at  the  boundary  of  this  region,  such  a  trajectory 
is  a  "barrier,"  across  which  the  time-to-go  is  discontinuous,  and  along 
which  the  gradient  is  infinite.  These  parameters  satisfy  the  relation 
P2  +  72  <  1  ,  the  equality  being  shown  as  ^  in  Fig.  2.4.  If  the 
parameters  fall  above  this  curve,  a  tangential  trajectory  cannot  arrive 
at  the  boundary  of  the  usable  part.  It  has  also  been  shown  that  if  0 
is  sufficiently  small,  the  capture  region  is  a  finite  curvilinear  tri¬ 
angle.  If  E  is  initially  outside  of  this  region,  he  can  avoid  capture 

indefinitely,  using  the  sidestep  tactic  of  the  bullfighter.  This_finite 

2  2  r  ’ 


illUCiilUVViJ  ,  — - -  — ^  )'■  „ 

capture  region  exists  (for  f?  +  72  <  l  )  U  0  <  7  sin  7  +  Vl-T  -  1  • 

*By  contrast,  we  have  seen  in  Sec.  2.2  that  the  y-axis  above  the  EL  is 
a  UL  ,  since  there  S+  =  S~  =  0  . 
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The  latter  function  appears  as  C  in  Fig.  2.4.  It  has  also  been  shown 

A 

that  the  point  A  occurs  on  the  barrier  (as  assumed  in  Ref.  1)  only  if 
c(p/2y)  +  7  >  0  ;  C  is  the  corresponding  locus  in  Fig.  2.4. 

These  three  loci  determine  two  regions  in  the  parameter  space,  but 
optimal  trajectories  which  fill  the  playing  space  have  been  found  in 
this  chapter  only  for  parameters  in  Region  I.  A  typical  set  of  trajec¬ 
tories  for  parameters  in  this  region  is  shown  in  Fig.  2.5(a).  Because 
the  time-to-go  is  discontinuous  (and  infinite)  at  the  barrier,  where  the 
point  A  is  located,  paths  "from"  A  must  instead  be  regarded  as 
starting  from  points  infinitesimally  close  to  A  .  From  this  range  of 
initial  points,  trajectories  depart  to  terminate  at  angles  in  the  half¬ 
open  interval,  0  <  Q  <  ©  .  We  follow  the  convention  that  when  the 


FIGURE  2.4.  Regions  in  Parameter  Space 
for  the  Classical  Game 
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(a)  Barriers  Closed 


FIGURE  2.5.  Trajectories  and  Isochrones 
in  the  Classical  Game 
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barr  rs  close,  E  follows  the  semipermeable  barrier  when  possible,  E 
preferring  an  early  safe  contact  at  6[ip  to  a  later  termination  at 

0=0.  P  ,  of  course,  turns  toward  E  whatever  E's  strategy,  be¬ 
cause  in  this  case,  S  has  the  sign  of  x  ,  and  P  will  reasonably 

prefer  a  safe  contact  to  none  at  all.  Thus,  in  this  case  we  see  that  B 

is  a  locus  of  infinite  discontinuity  in  the  time-to-go,  and  it  is  an 
optimal  path  which  ends  tangentially  at  Q 

up 

Contours  of  constant  V  can  be  easily  calculated  as  the  dashed 
lines  in  the  figure.  These  isochrones  are  tangent  to  the  barrier  only 
below  the  point  A  .  The  discontinuity  in  E's  strategy  for  initial 
conditions  near  S  at  the  point  A  is  the  cause  of  this  change  in  the 
isochrones . 

The  more  complex  diagram  in  Fig.  2.5(b)  applies  for  parameters  in 
Region  II  to  the  right  of  Cg  ,  and  this  configuration  is  essentially 
as  given  in  Ref.  1.  Because  the  strategies  are  not  known  at  all  points, 
however,  it  may  happen  that  some  of  the  trajectories  so  far  found  are 
not  optimal.  In  fact,  this  actually  occurs  for  a  small  range  of  para¬ 
meters,  as  will  be  shown  in  Chapter  III.  For  present  purposes,  we  may- 
say  that  the  trajectories  of  Fig.  2.5(b),  though  incomplete,  are  quali¬ 
tatively  representative,  for  all  parameters  right  of  Cg  and  above  C 
in  Region  II.  Of  course,  as  (0,7)  change  in  Region  II,  the  point  A 
will  be  displaced  along  the  barrier  and  the  barrier's  relative  size  will 
change,  but  this  does  not  affect  the  type  of  solution.  We  will  find, 
however,  that  the  reverse-chase  is  absent  for  parameters  close  to  the 
intersection  of  C  and  C  in  Region  II. 

The  isochrones  of  Region  II  are  more  complex  and  interesting  than 
those  of  Region  I.  The  barrier  appears  as  a  discontinuity  in  the  time- 
to-go,  but  this  discontinuity  decreases  to  zero  at  the  end  B  of  the 
barrier.  The  EL  ,  which  extends  the  barrier  to  the  y-axis  or  to  the 
capture  circle,  is  a  sharp  minimum  in  V  ,  in  the  sense  that  the  time- 
to-go  increases  linearly  in  both  directions  normal  to  the  EL  .  The 
normal  to  the  local  isochrone  is  thus  indefinite  on  the  EL  ,  reflecting 
the  fact  that  E  can  choose  between  two  strategies  on  the  EL 
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Similarly  we  have  found  that  when  the  barriers  are  open,  the  y-axis 
may  be  a  UL  or  a  PDL  ,  depending  on  whether  or  not  P's  switch  func¬ 
tion  is  continuous  across  it.  For  the  discontinuous  (PDL)  case,  the 
normal  to  the  isochrone  is  again  double-valued,  because  the  y-axis  is 
here  a  sharp  maximum  in  the  time-to-go,  V  . 

For  initial  conditions  on  the  open  barrier,  P's  strategy  is  a  turn 
away  (<p  =  -1)  unless  the  parameters 
these  parameters,  we  will  find  that  a 
for  which  Cp  =  +1  is  optimal  on  both 
is  an  optimal  trajectory  only  in  this 


are  below  C,  in  Fig.  2.4.  For 
4 

range  of  positions  on  S  exists, 
sides  of  the  barrier.  Thus,  S 
case,  and  when  the  barriers  close 
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Chapter  III 

2  2 

EXTENSIONS  TO  THE  CLASSICAL  RESULTS.  g  +  7  —  1 

The  parameter  inequality  in  the  chapter  title  holds  whenever  the 
barrier  exists  at  6up  •  But  when  the  right  and  left  barriers  do  not 
intersect,  strategies  for  x  in  the  shaded  area  of  Fig.  2.5(b)  must  be 
found.  For  parameters  above  C4  (see  Fig.  2.4),  including  those  to  the 
left  of  C3  ,  this  extension  is  not  difficult,  though  it  entails  treat¬ 
ment  of  "safe-contact”  strategies.  For  parameters  between  C2  and  C4 
in  Region  II,  however,  two  new  exceptional  lines  are  required  to  complete 
the  solution.  We  begin  by  considering  those  trajectories  which  touch  the 
capture  circle  at  the  lowest  corner  of  the  shaded  area. 


3.1  One-Dimensional  Motion 

We  have  seen  in  Fig.  2.5(b)  that  turn-away  trajectories  can  arrive 
at  the  negative  y-axis.  The  particular  path  which  contacts  the  capture 
circle  at  the  right  edge  of  the  shaded  region  is  specified  in  terms  of 
the  intercept  yx  and  the  time  from  tangency  to  this  intercept.  The 
equations  giving  these  two  quantities  are  r  =  P  ,  r  =  0  ,  which  re¬ 
quire  the  use  of  the  retrograde  trajectories 

x  =  -1  +  ct  -  (y1-n)st 

(1) 

y  =  st  +  (y1-7t)ct  . 


The  tangency  condition  is  then  expressed  by  the  equations 
i  -  ct  +  (y1-7t)st  +  4[ (y,-7i)2  -  P2]  =  0 

(2) 

(1-7)st  +  (y1-7i)  (ct-7)  =  0  • 

Eliminating  yx  gives  a  cubic  in  ct  ,  which  turns  out  to  be  the  square 
of 

v,1+t2-27ct  (1-ct)  -  (3(ct-7)  =  0  • 

This  expression  yields  ^  .  the  time  from  tangency  to  the  negative  y- 
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axis,  in  terms  ol  which  the  intercept  is 


yl  =  7T1  +  7-CT, 


<  -p  • 


The  angle  at  which  tangency  occurs  in  given  by  (1)  as 

ycTl-i 


tan  em  =  yCy,.^) 


(4) 


(5) 


7ST, 


The  supplement  of  this  angle  is  greater  than  0up  ,  as  can  be  shown 

by  using  (3)  and  (4),  together  with  the  definition,  c0up  =  7  •  That  is, 

E  is  physically  capable  of  maintaining  safe-contact  for  angles  greater 

than  d  but  such  a  strategy  is  not  optimal, 
m 

We  consider  next  the  conditions  at  tangency  with  the  aim  of  extend¬ 
ing  the  trajectories  backwards  from  this  point.  The  main  equation  in 
polar  coordinates  is 


min  max[V  r+V  0]  =  -1  . 

cp  lit  r  ° 


(6) 


where  the  polar  equations  of  motion  are  given  in  (11)  of  Sec.  2.3, 
r  =  -ce  +  ycty-e) 
e  =  -<p  +  i  [se  +  ysOl'-e)]  • 

i  _n  relates  t  to  0  (ambiguously),  while 
The  tangency  condition  r  -  0  relate  v 

(6)  reduces  to 


min  Vf-cp  +  |  [sS  +  7s(*-e)]}  =  -1 
cp  o  P 


(7) 


(8) 


where  the  angular  adjoint  is  identical  with  P's  switch  function; 

ve  =  vxy  -  Vyx  *  S  . 

The  minimizing  control  of  P  is  thus  q>  =  sgn  S  =  -1  ,  and  E's 
strategy  must  maximize  V  s  ((r-fl)  while  satisfying  f  =  0.  Since 
V  <  0  at  the  point  on  the  capture  circle  where  safe  contact  ends,  the 

e 


ambiguity  in  E's  control  is  resolved.  Thus, 


=  g  -  cos  1  (ce/7) 


(9) 
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is  E's  safe-contact  strategy  when  cp  =  -1  .  In  one-dimensional  (tan¬ 
gential)  motion,  we  find  that  E's  angular  rate  is 


e  =  i  +  (se  -  7 \/i-c2e/72)/g  ,  do) 

and  consequently, 


=  -1/e  =  - 


P+S0-7  V/1-c20/72 


(11) 


We  may  now  consider  the  trajectories  which  arrive  at  the  capture 

circle,  there  to  begin  a  safe-contact  path  which  departs  tangentially  at 

the  angle  0  given  in  (5).  Just  prior  to  contact,  the  main  equation 
in 

for  cp  =  -1  is  written 


max[V 

t  r 


r  +  vgn  =  -i . 


(12) 


Since  r  =  p  is  the  boundary  of  the  regular  region,  the  angular  gradient 
on  arrival  is  continuous.  The  maximization  of  (12)  implies  that 


c(i|r  -0) 


V 

r 


v/v^(V0/fJ)2 


(13) 


and  therefore,  using  (11)  for  V  /p  ,  (12)  can  be  solved  for 

6 

V  =  - SiZl - _____  .  (14) 

r  Vi-^e/r2  (p+s0-7  v/i-c20/72) 

This  finally  permits  (13)  to  be  solved  for  if  =  0-cos  1( cS /y)  =  i|/+  >  and 
therefore  r  =  r+  =  0  ,  so  that  trajectories  and  strategies  are  smooth 
on  arrival  at  the  capture  circle.  In  fact,  the  continuity  of  if  could 
have  been  anticipated  from  the  fact  that 

max(V  r  +V.0  )  =  max  V  0+  =  -1  , 
ijr“  r  9  i| 0 


i.e.,  because  E's  vectogram  is  circular,  H  has  a  unique  maximum 
relative  to  if  ,  such  that  if”  =  if+  .  A  typical  safe-contact  trajectory 
for  cp  =  -1  is  shown  in  Fig.  3.1. 
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3.2  Switch  Line  Replaces  Equivocal  Line 


For  parameters  slightly  to  the  left  of  C3  in  Fig.  2.4,  the  point 
A  falls  between  the  end  B  of  the  barrier  and  the  negative  y-axis. 

The  points  B  and  A  are  connected  by  a  switch  line  (SL)  ,  on  which 
S  =  0  ,  and  A  is  connected  to  the  capture  circle  or  to  the  y-axis  by 
an  EL  , 


The  disappearance  of  the  EL  occurs  when  the  point  A  falls  on  the 
capture  circle.  This  is  easily  derived  as  the  locus 

c(g/7)  =  l-p2/2  ,  (1) 

which  will  be  labelled  in  the  parameter  space  drawing  at  the  end  of 

this  chapter  (Fig.  3.9).  For  parameters  to  the  left  of  this  line,  retro¬ 
gressive  trajectories  to  A  are  interrupted  by  the  capture  circle,  and 
the  point  B  is  connected  to  the  capture  circle  by  the  SL  , 


The  retrogressive  switch  time  is  determined  as  a  function  of  arrival 

angle  0  : 

o 


S 


V 

x 


y-v  x 

y 


ce-cie  +t) 

o  o 

ceo-7 


=  0 


or 


t  =  2(n-0  ) 
s  o 


(2) 


The  arrival  angle  91  ,  as  defined  in  Fig.  3.2,  is  given  by  equating 

T  =  t  =  B/7  ,  or  9.  =  it-p/27  .  For  all  arrival  angles  between  0  and 
s  A  1 

0^  ,  A  serves  as  a  dispersal  point  for  E  ,  and  it  is  interesting  that 
in  this  case  trajectories  can  arrive  at  A  from  two  directions. 

At  the  end  of  the  barrier,  the  trajectory  has  the  slope 

l-x-7c  (0  +t) 

_  _ 2  (3) 

dx  ”  y-7s(eo+tT”  ’ 

where  the  point  B  is  evaluated  using  C0o  =  7  together  with  (1)  to 
give,  as  in  (14)  of  Sec.  2.3, 
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T 


x 

y 


(4) 


Si 

dx 


/1-72 


This  implies  that  dr/dt  =  0  at  this  point.  The  slope  of  the  SL 
is 


Si 

dx 


S 


y  +y  (dt  /d0  ) 

9  t  s  o 

x  +x  (dt  /d0  )  ’ 
6  T  S  o 


(5) 


where  dr  /d 9  =  -2  ,  according  to  (2),  and  the  partial  derivatives  fol 

s  o 

low  immediately  from  the  retrogressive  equations  for  the  trajectories. 
Substitution  and  simplification  yield  the  result 


Si  ^-T2  (6) 

dx  s  "  7  ' 


which  with  (4)  proves  that  the  SL  is  tangent  to  the  trajectory  at  this 
point. 

At  the  point  A  ,  however,  where  0/7  =  2(n-91)  ,  the  trajectory 
slope  is  found  in  term3  of  the  arrival  angle  0j  , 


d* 

dx 


T 


c291-7c61 
-s291+7901  ' 


(7) 


This  is  the  same  result  given  for  the  slope  of  the  SL  at  this  point, 
as  can  easily  be  verified.  The  second  extreme  trajectory  from  A  ar¬ 
rives  at  the  capture  circle  at  9=0,  and  its  slope  at  A  is 

dy  _  c(3/7)  _  C26l  (8) 

dx  “  s (g/7)  ”  -s291 

The  EL  slope  at  the  point  A  requires  use  of  the  intermediate 
in  the  equations  of  motion,  and  the  result  is  expressible  as  a  function 
of  7  and  0  .  It  has  not  been  possible  to  derive  analytical  conclu¬ 

sions  from  the  expression.  However,  numerical  studies  for  a  range  of 
parameters  left  of  Cg  in  Region  X  show  that  the  EL  slope  is  not  equal 
to  either  of  the  values  (7)  or  (8). 

The  simplest  way  of  expressing  the  SL  is  as  a  function  of  arrival 
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angle  0q  ,  by  using  the  value  x  =  2(n-$°)  in  the  retrogressive  equa¬ 
tions;  i.e., 

x  =  ae  [286  -  p  +  2?(n-9  )] 

o  o  O 

r  (9 

y  =  -ce  [2S3  -  p  +  2  y  ( jt-e  )] 

O  o  O 

where  0  <  6q  <  cos_1y  (Cf.  (14)  of  Sec.  2.3). 

The  qualitative  appearance  of  the  trajectories  for  parameters  left 
of  C5  ,  which  is  given  by  (1),  is  as  shown  in  Fig.  3.3.  Since  S 
passes  smoothly  through  zero  for  trajectories  intersecting  the  SL  ,  it 
follows  that  E  is  running  toward  P  when  the  SL  is  crossed.  Ex¬ 
cept  for  the  absence  of  an  EL  ,  this  configuration  is  very  much  like 
that  of  Fig.  3.2.  We  note  that  in  both  of  these  figures,  a  local  maxi¬ 
mum  in  the  optimal  time-to-go  occurs  at  a  point  on  the  underside  of  the 
barrier  roughly  one-fourth  of  the  way  to  the  end  of  it.  The  loci  of 
constant  V  can  be  regarded  as  elevation  lines  on  a  contour  map  of  the 
playing  space.  This  concept  will  be  helpful  in  understanding  some  re¬ 
sults  to  be  found  in  the  next  section. 


3'3  Switch  Envelope  and  Pursuer's  Dispersal  Line 

For  parameters  in  the  lower  portion  of  Region  II,  it  is  found  that 
IB  is  not  necessarily  a  locus  of  discontinuity  in  P's  strategy.  That 
is,  over  a  certain  interval  of  IB  ,  the  switch  function  is  positive  on 
both  sides  of  it,  such  that  when  E  is  under  the  barrier  here,  P 
lunges  toward  him.  The  verb  is  appropriate  because  this  is  always  a 
very  brief  maneuver,  soon  followed  by  a  turn-away  and  subsequent  safe- 
contact.  As  mentioned  in  Sec.  2.5,  when  cp  changes  sign  across  B  , 
the  barrier  is  not  an  optimal  path  but  is  only  a  locus  of  initial  condi¬ 
tions  for  the  turn-away  maneuver.  But,  when  parameters  are  such  that 
the  lunge  region  exists,  we  can  expect  P’s  switch  line  to  lie  to  the 
right  of  IB  (which  now  i^  an  optimal  path),  somewhat  as  shown  in  Fig. 
3.4.  We  consider  a  comparative  test  of  the  three  retrograde  paths 
shown  in  this  figure. 
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At  the  point  ,  when  the  lowest  path  reaches  the  barrier, 

it  is  found  that  if  P  were  to  switch  to  q>”  =  +1  ,  thereby  forcing  E 
to  switch  to  a  strategy  i|/  perpendicular  to  £  ,  then  for  x  on  S  , 

W+.x~  >  -1  .  (!) 

The  inequality  Implies  that  such  a  switch  is  to  E's  advantage  and 
hence  that  <p  =  +1  is  not  optimal  for  such  points.  It  may  also  happen, 
on  the  other  hand,  that  at  G  on  a  second  path  T„  ,  this  double 
switch  in  strategies  leads  to  the  inequality 

W+-x  <  -1  ,  (2) 

which  means  that  P  can  force  E  to  cross  the  'V*  contours  faster  than 
the  optimal  rate.  Therefore,  these  contours  are  not  correct  at  this 
point  and  E's  strategy  cannot  be  given  by  W+  in  this  region,  where 
P  lunges  toward  E  . 

If  these  inequalities  hold,  then  at  some  intermediate  point  G  of 
the  barrier,  we  must  find 

W+.x"  =  V*x"+V*y"  =  -1  ,  (3) 

where  the  left  side  can  be  expressed  as  a  function  of  0  ,  by  using 

x+  =  x  to  provide  tb  and  tq  in  terms  of  0q  .  Specifically,  we 
have 


x  =  1  -  ex.  +  (3-71.  )s(0  +T.) 

b  b  up  b 

y  =  St.  +  (3-7t.  )c(S  +t  ) 
b  b  up  b 

X+  =  -1  +  CT  +  9s  (0  -T  )  -  7T  sCllt  -T  ) 
O  O  O  0  o  o 

y+  =  sto  +  3c(0o-to)  -  7toc(Vo-to)  , 


(4) 


(5) 


where  \|r  =  0  -  cos  1  (c0  /y)  and  6  =  cos_17  .  Finally,  the  ad- 

o  o  o  up 

joint  components  in  (3)  are  easily  found  in  terms  of  \|/o  and  the  retro¬ 
gressive  time  from  tangency,  tq  ; 

V+  =  s(f  -t  )/D 

X  o  o 

+  <6> 
Vy  =  C(VT0)/D 
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D  =  V1-c20o/t2 (3+s0o~7  v 1-c 29o//J )  . 

By  numerical  search  methods,  we  determine  the  parameters  0q  ,  tq 

and  x  which  specify  the  endpoint  G  of  the  SE  .  Of  course,  when  the 
b 

condition  (3)  does  not  hold  for  any  point  on  55  ,  the  parameters  (3,7) 
are  above  the  locus  C4  of  Fig.  2.4.  Since,  whenever  (3)  holds,  the 
barrier  always  has  the  qualitative  form  of  Fig.  3.4,  the  point  G  is 
the  lowest  point  of  a  "lunge"  region.  We  then  must  find  the  switch  line 
which  begins  at  G  and  which  separates  this  lunge  region  (9  =  +1) 
from  the  turn-away  region  (9  =  -1)  . 

To  develop  generating  equations  for  this  line,  we  first  consider  a 
group  of  trajectories  parallel  to,  and  right  of,  55  .  At  the  point  G  , 


FIGURE  3.4.  Detail  Near  Barrier  When  SE  Exists 
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we  have  (3)  together  with  the  main  equation  on  the  prior  side,  W  •£  - 

-1  ,  which  combine  to  give  the  slope  of  the  switch  line  as  a  function  of 


position, 


V  -V 

d£  _  _  _x _ x 

dx  "  v_-V+ 


(7) 


On  3  ,  of  course,  |W|  -* »  ,  and  so  the  slope  of  the  switch  line  is 
the  same  as  the  slope  of  the  local  trajectory,  3  .  As  the  SL  is  fol¬ 
lowed  away  from  3  ,  the  incoming  trajectories  must  continue  to  be  tan¬ 
gential  to  it.  This  is  so  because,  if  they  encountered  this  switch  line 
at  a  non-zero  angle,  and  if  P  chose  to  delay  his  turn-away,  E  would 
"chatter"  between  his  strategies  f  and  *+  ,  thereby  remaining  on  the 
line.  Since  E’s  vectogram  is  a  circle,  however,  this  chatter  strategy 
would  actually  reduce  his  speed  and  hence  would  not  yield  dV/dt  =  -1  , 
as  did  9  =  +1  and  *r  «  *“  .  Hence,  the  switch  line  is  the  envelope  of 
the  prior  paths,  and  consequently  it  is  called  a  "switch  envelope"  (BE). 


While  they  are  similar  in  some  superficial  respects,  the  SE  dif¬ 
fers  in  an  important  way  from  the  EL  described  in  Chapter  II.  The 
switch  function  goes  smoothly  to  zero  at  the  endpoint  B  of  3  ,  which 
is  then  tangentially  extended  by  the  EL  .  It  was  found  that  S  =  0 
on  the  incoming  paths  to  the  EL  ,  while  S+  >  0  on  the  departing  tra¬ 
jectories.  This  meant  that  E  could  choose  to  kee£  S  =  0  by  ^running 
radially,  which  gave  P’s  control  as  the  intermediate  value,  $  .  By 
contrast,  neither  s“  nor  S+  is  zero  along  the  SE  ,  but  it  is  known 
that  the  component  of  W  along  the  SE  is  continuous,  so  that  (3) 
holds  along  it.  This  implies  that  P’s  control  along  the  SE  is 

=  +i  ,  and  hence  E’s  arrival  control  is  given  by  (3)  in  terms  of 

e  and  T  •  The  expression  (7)  for  the  SE  slope  is  then  solved  for 

O  o 

the  derivative  d6o/dt0  • 


As  described  in  detail  in  Appendix  A,  the  retrogressive  integration 
of  the  SE  is  valid  only  so  long  as  the  curvature  of  the  SE  exceeds 
that  of  the  incoming  paths.  This  curvature  condition  is  equivalent 
condition  of  derivatives  of  f  .  On  the  incoming  paths,  *  "  “l  • 
while  on  the  SE  ,  if  is  expressed  as  a  function  of  0Q  and  fQ 
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and  its  time  rate  of  change  is  given  in  terms  of  these  coordinates.  In 
other  words,  <jr  is  discontinuous  at  tangency,  except  at  the  upper  end 
of  the  SE  . 

At  the  coordinate  Qq  ,  where  these  curvatures  are  equal,  the  SE 
is  tangentially  extended  by  a  PDL  .  That  this  should  be  so  is  suggested 
by  Fig.  3.5,  where  the  cotangency  of  PDL  ,  SE  and  trajectory  at  their 
meeting  point  is  implied  by  a  finite  curvature  condition  on  the  trajec¬ 
tory  at  the  junction. 


FIGURE  3.5.  Detail  Near  Junction  of  SE  and  PDL 

The  PDL  is  simpler  to  construct  than  the  SE  ,  though  the  PDL 
is  given  in  terms  of  SE  coordinates,  which  must  therefore  be  stored  as 
they  are  computed.  The  retrogressive  trajectories  of  either  side  must 
yield  the  same  point  on  the  PDL  ,  and  the  time-to-go  must  be  the  same 
on  the  two  paths.  Denoting  the  side  near  the  barrier  by  the  super¬ 
script  "+"  (corresponding  to  the  lunge  strategy,  o  =  +1),  we  have 
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(8) 


CiX+teo'V  =  dx”(91-'t1> 

dy+(9o,T2)  = 

dV+(9l’T2,0o)  =  dV_(T1>  * 

where  the  various  angles  and  times  are  as  defined  in  Fig.  3.5.  The  ap¬ 
propriate  differential  equations  required  in  generating  the  PDL  are 

symbolically  expressed  using  (8),  with  Q  taken  as  a  convenient  inde- 

o 

pendent  variable  (subscripts  denote  partial  differentiation): 


+ 


+  x 


T  d  Q 

o 


e  d  e 


dx. 


t  de 


ye  + 


-  d9i  -  dTi 

ye  de  +  yT  de 

o  o 


(9) 


At  given  values  of  the  dependent  variables,  the  gradients  W+  and  W 
are  known  for  use  in  the  slope  equation,  (7).  This  allows  computation 
of 


dx2  yg-Xg(dy/dx) 

dtP  =  “  +  ..  ,  > 

0  y^-x^tdy/dx) 


(10) 


and  a  straightforward  inversion  of  (9)  for  the  remaining  two  derivatives. 

Initial  values  at  the  point  D  are  t  =  0  ,  0=0  and  t.  =  T  , 

2  1  o  1  o 

the  SE  integration  having  provided  0  and  t 

o  o 

The  upper  end  of  the  PDL  occurs  at  the  barrier,  and  these  curves 
Join  here  tangentially  owing  to  the  infinite  magnitude  of  the  gradient 
vector  vfy+  at  the  barrier.*  That  is,  the  PDL  slope,  as  given  by  (7), 
is 


51  =  _ 

dx 


v“-v+ 

X  X 

v”-v+ 
y  y 


(li) 


*For  some  parameters  near  the  lower  edge  of  Region  II,  the  PDL  actually 
extends  outside  the  shaded  area  of  Fig.  2.5(b).  This  illustrates  that 
unless  the  entire  space  is  filled  with  trajectories,  none  of  them  can 
strictly  be  called  "solutions." 


42 


which  is  equal  to  the  slope  of  the  adjacent  barrier,  as  obtained  using 
(2)  of  Sec.  2.3, 


*1  =  _ 

dx 


(12) 


Since  3  is  a  trajectory  on  the  "+"  side  of  the  PDL  ,  the  optimal  time- 
to-go  is  a  local  maximum  somewhere  to  the  right  of  the  tangential  junc¬ 
tion  of  PDL  and  3  ,  as  will  be  shown  in  Fig.  3.7. 


3.4,  Focal  Line 


For  parameters  which  are  near  the  right  corner  of  Region  II,  numeri¬ 
cal  generation  of  the  SE  shows  that  at  some  point  below  where  the  cur¬ 
vature  condition  denoting  the  start  of  the  PDL  is  satisfied,  the  SE 
derivative  vanishes,  or  =  0  .  At  this  point  on  the  SE  ,  the 

velocities  obtained  by  two  sets  of  controls  are  identical;  i.e.,  at  a 
single  point,  incoming  and  outgoing  velocities  are  equal,  x(cp  ,\|r  )  = 
x(<p+,i|r+)  .  The  point  where  this  occurs  is  found  to  be  the  lower  end  of 
a  "focal  line"  (FL) ,  along  which  E  runs  regardless  of  P's  strategy. 
Furthermore,  the  FL  is  an  envelope  of  arrival  paths  on  both  sides. 

The  incoming  paths  are  tangential  because,  if  P  maintains  his  arrival 
strategy,  E  must  do  so  as  well,  in  order  to  keep  dV/dt  =  -1  .  As  in 
the  case  of  the  SE  ,  both  strategies  are  continuous  when  a  path  joins 
the  FL  ,  and  therefore  the  paths  themselves  have  continuous  slopes  and 
are  tangential  at  the  FL  . 

The  FL  is,  like  the  SE  ,  a  local  (sharp)  maximum  or  crest  in  the 
time-to-go  contours,  and  for  an£  choice  of  cp  ,  E  can  stay  on  this 
crest.  Progress  down  the  line  will,  however,  be  maximized  if  P  chooses 
an  extreme  control,  cp  =  il  .  As  indicated  in  Fig.  3.6,  the  FL  condi¬ 
tion  yields  two  equations  in  t|c  and  ,  which  are  E's  controls 

corresponding  to  cp  =  -1  and  +1  ,  respectively.  Using  the  equations, 
then, 


x  =  -y  +  7S'I'+  =  y  + 

y  =  -1  +  x  +  7c4c+  =  -1  -  X  +  7C\|r  . 


(1) 
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(x,y> 


FIGURE  3.6.  Focal  Line  Strategies 


These  give  E's  strategy  as  a  function  of  both  E's  position  on  the 
FL  and  P's  strategy  (cp  =  +1)  ,  as  follows: 

ysty  =  ipy  +  (x/r)  v'T? 

<2> 

7Cf  =  -ox  +  (y/r)  V  TT-ir  . 

Using  these  expressions  in  (1)  gives  the  retrogressive  differential  equa¬ 
tions  for  the  FL  as  the  simple  pair. 


#  =  -(x/r)  v  ^-r2 
?  =  1  -  (y/r)  v  T2-!-2  . 


(3) 


Neither  strategy  is  explicit  in  (3),  but  =  il  and  \|f(tp)  Is  given  by 

(2). 
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It  is  necessary  to  have  the  adjoints  along  the  FL  as  well,  since 
they  are  required  in  the  subsequent  integration  of  the  SE  .  We  know 
sV/Vx  =  c\|//V^  on  both  sides,  for  either  choice  of  y  ,  and  the  main 
equation  gives 

V^x  +  V^y  =  (xs\jc+yc\|/)/D  =  -1  . 

This  can  be  solved  for 


D  =  -^[(y/r)  V  ^-r2  -  y2  +  r2  -  cpx] 


and  this  denominator  (the  inverse  of  the  gradient  magnitude)  is  positive 
for  x  on  the  FL  .  The  direction  of  the  gradient  vector  on  each  side 
is  parallel  to  E's  strategy,  (2). 


The  FL  has  features  unlike  those  of  any  other  type  of  exceptional 
curve  found  in  this  game.  The  equations  (1)  show  that,  to  traverse  the 
FL  ,  P  can  choose  either  extreme  rate  of  turn,  and  E’s  strategy  de¬ 
pends  on  this  choice,  according  to  (2).  Further,  P  can  switch  strate¬ 
gies  ad  libitum  between  +1  and  -1  ,  all  with  no  change  in  the  rela¬ 
tive  trajectory,  provided  that  E  reacts  immediately.  Examination  of 
the  vector  velocities  in  Fig.  3.6  shows  that  no  intermediate  turn  rate 
can  produce  the  same  relative  motion  as  cp  =  il  .  These  comments  are 
meant  to  emphasize  the  distinction  between  "chatter"  on  the  EL  ,  and 
"switching"  on  the  FL  .  A  further  observation  is  that  the  switch  func¬ 
tion  is  never  zero  on  the  FL  ,  while  it  remains  at  zero  for  motion 
along  the  EL  . 


The  upper  end  of  the  FL  occurs  when  its  curvature  is  the  same  as 
that  of  an  incoming  path.  This  curvature  criterion  is  simply  written  as 
a  rate  of  change  of  i|r  along  the  FL  ,  or 


d£ 

dx 


(4) 


If  this  equation  is  first  satisfied  for  cp  =  -1  ,  for  example,  the  in¬ 
coming  trajectories  to  the  right  have  the  same  curvature  as  the  FL  , 
and  hence,  as  shown  in  Fig.  3.7,  the  FL  is  extended  by  an  SE  .  When 
the  time  derivatives  are  written  as  in  (3),  we  find  the  criterion  as 
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.  2 

ry  -  v  r-r2,  (2r  -Kpx)  =  0 


(5) 


where  rp  =  +.1  .  For  the  many  cases  of  the  focal  line  which  have  been 
analyzed  numerically,  (5)  was  first  satisfied  for  <p  =  “1  ,  so  that  the 
FL  is  extended  by  a  SE  .  This  means  that  a  "reverse  SE  with  tra¬ 
jectories  tangent  on  the  right,  and  for  which  (5)  would  be  first  satis¬ 
fied  for  cp  =  +1  ,  is  apparently  never  realized. 

Extension  of  the  FL  by  an  SE  is  straightforward  except  for  the 
initial  step  in  the  integration.  When  =  0  ,  using  the  notation  of 

Fig.  3.8,  the  indeterminacy  in  the  slope  of  the  SE  is  evaluated  by  ex- 


FIGURE  3.8.  Detail  Near  Junction  of  SE 
pansion  in  series.  Thus,  the  slope  of  the  SE  is 


and 


FL 


d£  _  £2  =  i-x-ycy 


dx 


S 


y-7SV 


where 


:Of-VV  =  1  +  "  CVf 


-  r‘/7] 


(6) 


(7) 


i7 


The  expression  (7)  follows  from  the  SE  condition  and  the  trajec¬ 
tories,  which  here  are  given  as  functions  of  ^VT1^  • 


x  =  -1  +  (xf+l)cx1  -  yfst1  "  5't1s(\(ff-T1) 


(8) 


y  =  yfcxx  +  (xf+l)st1  -  7T1c('l<f-T1)  • 

The  subscript  f  denotes  a  functional  dependence  on  xf  .  Then  the 
slope  (6)  is  equated  to  the  bilinear  form, 

yT  +yT<dVdTf> 

dy  _  f _ 

dx  “  x  +x  (dx,/dx_)  ’ 

T,  T  1  f 


which  is  inverted  to  yield 

Hi  _  _ 


x  f  -y  £  v 

TV  T.  _  A(x.  ,  T_) 

I _ X  &  l  I 

p  “  o-  B  (x..  ,  T  -  ) 

x  y  -y  x  If 


(9) 


It  is  easy  to  see  that  A  =  B  =  0  when  =  0  ,  since  then  all  veloci¬ 
ties  are  equal;  e.g.f  x^  =  ~  ^  •  With  close  attention  to  detail) 

it  can  also  be  shown  that  the  first  partial  derivatives  of  A  and  B 
vanish  at  this  time.*  Consequently,  (9)  must  be  expressed  as 


2  2 
A  dx.,  +2 A  ^  dx.  dx„+A  dx„ 

dTl  Vl  1  Vf  1  f  TfTf  f 

dTf  B  dx^+2B  dt  dT  +B  dxj 

Yi  if  f f 


(10) 


It  can  then  be  shown  that,  at  x^  =  0  ,  ATlTf  =  AXfXf  =  BXjXf  ~  0  > 
so  that  (10)  becomes  linear;  the  solution  is 


Hi 

dt* 


A  +2B 

Vl  Vf 


Vl 


(11) 


v° 


where 


*The  upper  end  of  the  FL  is  implied  by  (5)  of  Sec.  3.4,  and  this  func¬ 
tion  is  a  factor  of  dB/dx^  . 
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t  =  (A^-A^) (1-yv)  +  (B1~B2)xv 

B  -  (A  -A  )  (1-yv)  +  (B  -B  )xv 
T,  T  O  2  0  2 


2B  =  (r2-y2+x+yv)  [" -^(2y-r2v)  +  4“  (v  "  “2"2  '  *2  ) 
ll  f  Lr  r  v  \  r  v  r  /J 


Ao  =  2  XV  ‘  y 


l  .  2  2  4.  ,  4 

A1  =  y  -  Jxv  +  (y  x  -r  )/r  v 


A2  =  J(xv+y)C2  -  (yv-x)C2  +  Jy  -  xv 


„  3  , 

b  =  «  yv  +  x  -  i 

O  £ 

2  4 

B  =  -x  +  J(l-yv)  +  y  xy/r  v 

B„  =  i(yv-x)C2  +  (xv+y)C  +  i(l-x)  +  yv 

2  i  £ 

2  2 

C1  =  1  +  x/r  -  y/r  v 

2  2  2  2  6  3 

C2  =  (y+x/v)/2r  -  (7  -2r  >(xv-y)  /4r  v 

v  mJl y2/eb-l  , 

and  x,y  and  r  =  v^+y2  denote  values  at  the  junction  of  FL  and  SE 

The  SE  continues  to  extend  upward  until  a  curvature  condition 
marks  its  tangency  with  a  PDL  ,  as  in  Sec.  3.3.  The  PDL  is  then  found 
to  continue  until  its  tangency  with  8  .  This  final  integration  ends  at 
the  point  where  8  is  a  trajectory  connecting  the  top  of  the  PDL  with 
the  bottom  of  the  SE  ,  as  shown  in  Fig.  3.7.  Also  shown  in  that  figure 
are  the  contours  of  constant  V  ,  and  as  mentioned  earlier,  a  local 
maximum  in  time-to-go  must  exist  at  a  point  on  8  which  is  to  the  right 
of  the  PDL  junction. 
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Parameter  Space;  Regions  I  and  II 

The  current  version  of  the  parameter-space  diagram  is  shown  in  Fig. 

3.9.  The  locus  C.  is  the  line  above  which  the  SE  and  the  attendant 
4 

PDL  are  not  present  beneath  the  barrier.  It  is  determined  numerically, 
by  seeking  to  satisfy  (3)  of  Sec.  3.3  at  a  fixed  7  ,  and  for  variable 
(3  .  The  curve  C.  represents  parameters  for  which  the  dispersal  point 
A  occurs  on  the  capture  circle,  and  is  easily  found  as  (1)  of  Sec. 

3.2.  Finally,  the  small  triangular  zone  bounded  by  and  Cg  in¬ 

cludes  those  parameters  for  which  the  SE  is  interrupted  by  the  FL  ; 

the  locus  0  is  also  determined  numerically, 
b 


In  completing  Region  II  of  the  parameter  space,  diagrams  of  optimal 
paths  and  isochrones  have  been  given  in  Figs.  3.2  and  3.3,  for  para¬ 
meters  left  of  C3  .  For  parameters  below  C 4  ,  the  gross  features  of 
the  trajectories  in  the  area  under  the  barrier,  including  the  SE  ,  FL 
and  PDL,  have  been  sketched  in  Figs.  3.4  through  3.8. 
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We  note  again  that  if  the  barriers  do  not  intersect  on  the  y-axis, 
the  barrier  is  not  an  optimal  path  unless  the  SE  exists.  In  this  case 
(8,7  between  C  and  C  ),  only  that  segment  of  %  between  the  SE 
and  the  PDL  is  an  optimal  path.  The  remainder  of  it  is  a  locus  (de¬ 
scribed  with  cp  =  +1  )  of  initial  conditions  for  which  P's  optimal 
strategy  is  a  turn-away  (cp  =  —  1  >  . 
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Chapter  IV 


FURTHER  EXTENSIONS,  p2  +  y'  >  1 

The  classical  results  of  the  homicidal  chauffeur  game  have  been 
extended  in  Chapter  III  to  all  parameters  satisfying  the  inequality 
t  <  1  !  and  we  have  further  subdivided  this  area  into  Regions  I 
and  II,  according  to  whether  or  not  the  capture  region  is  finite.  The 
approach  taken  in  this  chapter  is  bnsieally  one  of  continuation  of  the 
loci  in  Fig.  3.9.  Beginning  the  extension  at  the  locus  C^  ,  we  will 
work  upward  in  parameter  space,  introducing  new  complexities  as  required. 
In  this  way,  previously  defined  loci  are  followed  to  their  ends.  Other 
loci  will  also  enter  the  discussion,  and  it  will  be  found  convenient  to 
separate  the  parameter  space  p2  +  >  1  into  three  distinct  regions. 

4 . 1  A  Mew  Barrier 

The  equations  describing  the  barrier  in  Region  I  were  derived  in 
Sec.  2.3,  starting  with  the  main  equation  in  the  terminal  payoff  form, 

min  max[v  x  +  V  y]  =  0  • 

<p  t  x  y 

2  2 

If  this  equation  is  considered  to  hold  at  r  =  (3  ,  but  with  P  +  y  >  i 

we  have  seen  in  (12)  of  Sec.  2.3  that  9  cannot  end  at  0up  with  tan¬ 
gential  velocity.  Neither  can  both  odjoints  be  zero.  The  remaining 
possibility,  if  this  equation  holds,  is  that  an  equilibrium  point  may 
exist  where  x  =  y  =  0  .  For  this  static  interpretation,  the  relative 
velocity  equations  are,  with  P  turning  toward  E  • 

pee  -  jsi/  =  o 

1  -  pse  -  7ci)r  =  0  . 

These  con  be  solved  for  the  equilibrium  position  and  E's  corresponding 
control, 


The  first  of  these  equations  implies  that  P  <  1  +  7  in  order  for  the 
equilibrium  to  exist.  But,  we  will  see  that  a  barrier  is  present  for 

>  l  only  when  |3  is  much  less  than  1  +  7  ,  so  that  the  right 
and  left  barriers  intersect  ahead  of  P  .  An  equilibrium  strategy  for 
P  (for  which  the  time-to-go  is  infinite)  can  be  minimizing  only  if  the 
capture  region  Is  bounded,  since  otherwise  P  could  certainly  choose  a 
better  strategy. 


When  the  retrogressive  equations  are  solved  with  the  appropriate 
boundary  conditions,  the  new  barrier  equations  are 


x  =  1  -  CT  +  pis  (t)  +t)  -  ns(\|/eti+-r) 

y  =  st  +  0c(Sgq+T)  -  7tc4eq+t)  . 


(3) 


The  parameters  for  which  the  right  and  left  barriers  intersect 
tangentially  on  the  y-axis  can  be  obtained  us  in  Chapter  II,  by  eliminat 
ing  T  between  the  equations  x((3,7,t)  =  0  ,  and  S(t3,7,x)  =  0  .  The 
result  is  p(7)  ,  the  extension  of  C 2  lying  to  the  right  of  , 
which  is  given  by  the  following  rather  unappetizing  implicit  function  of 
the  parameters: 


1  +  (5c(^-0eq)  -  s*eq  -  7<£  -  Veq>  =  0  • 

In  this  equation  0  and  if  are  first  quadrant  functions  of 

GCJ  GCJ 

defined  in  (2). 


(4) 

(P>?) 


The  switch  function  is 

S  =  7  -  c(tfi  +t)  , 

eq 

and  since  Ye  >  0eq  >  Sup  =  cos_17  ,  this  function  is  positive  for  all 
T  yielding  x  >0  in  (3).  That  is,  the  barriers  never  end  before 
meeting  on  the  y-axis  when  the  parameters  fall  below  the  locus  C2 
given  by  (4). 

The  trajectory  given  by  (3)  differs  in  three  respects  from  the 
barrier  trajectory  given  by  (7)  of  Sec.  2.3: 
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1)  Thu  arrival  velocity  is  zero  in  both  radial  and  tangential 
directions,  though  the  equilibrium  point  is  reached  from  any 
point  on  the  barrier  in  finite  time; 

ii)  The  new  barrier  is  not  tangential  to  the  capture  circle  at 
arrival; 

ili)  The  new  barrier  contacts  the  capture  circle  at  the  angle 

0  >  6 

eq  up 

This  new  barrier  is  the  outer  edge  of  an  angular  interval  of  the 
capture  circle  on  which  safe-contact  motion  is  optimal,  with  P  turning 
toward  E  .  Details  of  the  relative  motion,  including  the  influence  of 
an  evader's  "dispersal  point"  A  ,  are  given  in  the  next  section.  The 
normalized  gradient  vector,  which  is  parallel  to  E's  optimal  direction, 
is  normal  to  the  new  barrier,  as  was  the  case  in  Region  I.  The  gradient 
rotates  at  unit  rate  with  motion  of  x  down  the  barrier,  until  arrival 
at  the  equilibrium  point. 


4.2  Evader's  Dispersal  Point  and  Dispersal  Line 

When  the  parameters  are  outside  the  locus  ,  but  below  C ^  ,  the 

point  A  (and  its  image  across  the  y-axis)  is  Inside  rather  than  on  the 
new  barrier.  This  point  is  located  by  the  same  equations  found  in 
Chapter  I, 


xA  =  1  -  c(j}/7) 
yA  =  s((3/7)  , 


(1) 


and  trajectories  from  this  evader's  dispersal  point  are  expressed  as 

functions  of  the  arrival  angle  @  <  0  <0  and  the  retrogressive 

—  o  —  up 

time,  t  : 


x  =  1  -  ct  +  (0-7t)s (6+t) 
y  =  st  +  (p-7T)cte+t)  . 


(2) 


The  largest  angle,  .  is  given  by  solving  the  simultaneous  equations 

r((j  lT  )  =  (3  ,  ,  T,j,)  =  0  ,  as  detailed  in  Appendix  B-l,  and  as  shown 

in  Fig.  4.1. 
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This  figure  also  shows  that  A  is  the  near  end  of  an  evader's 
dispersal  line  (EDL)  extending  to  the  y-axis,  there  meeting  the  image 
EDL  from  the  other  disp.-rsal  point.  This  EDL  is  a  locus  of  starting 
points  for  which  E  must  choose  between  two  equally  optimal  strategies, 
labelled  "+"  and  for  reference  purposes  in  Fig.  4.1.  On  the 

side,  E  runs  along  the  tangent  to  P's  minimum  turn  circle,  and  the 
game  ends  as  a  straight-line  chase.  The  other  choice  of  strategy  from 


FIGURE  4.1.  EDL  from  Point  A  in  Region  III 


the  EDL  has  E 
ture  circle  near 
the  safe-contact 


running  straight  until  tangential  arrival  at  the  cap- 
the  equilibrium  point.  Subsequent  motion  here  is  of 
type  until  0  decreases  to  t?_  >  where  E  leaves  the 
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circle  lor  a  short  interval  beiore  termination  at  8q  •  It  Is  easily 
shown  that  P's  strategy  is  the  same  for  either  choice  of  E's  (i.e., 
S=Vy-Vx>0  if  x  >  0  ,  for  either  W  or  vV  ) .  The  -  tra¬ 
jectories  encounter  the  y-axis  as  in  Chapter  II,  with  P's  switch  func¬ 
tion  going  smoothly  to  zero  along  this  line.  On  the  barrier  side  of 
the  EDL  ,  though,  the  switch  function  is  positive  as  x  approaches  the 
brief  y-axis  segment  of  the  EDL  shown  in  Fig.  4.1.  While  it  can  be 
shown  only  numerically,  there  exists  a  discontinuity  in  both  W  and  S 
a.  toss  this  line  segment,  such  that  P's  strategy  depends  on  E’s 
choice. 

Referring  to  Fig.  B-l,  the  discontinuity  in  the  gradient  across  the 
EDL  gives  an  equation  for  the  local  slope  as 

V'-V+ 

d£  =  _  x  x  _  (3) 

dx  ”  v"-V+ 

y  y 


Next,  in  equating  x+  =  *“  *  differential  changes  lead  to  relations  be¬ 
tween  dd2,dt2,dy1,dT1  .  The  appropriate  independent  variable  is  yx  , 
since  the  other  variables  are  either  too  nearly  static  or  are  not  clearly 
monotonic.  Thus,  equating  differential  positions  on  either  side,  we 
derive  simultaneous  differential  equations, 


+  d02  +  dT2 

X  .  +  X  =  X  +  X 

0  dy  T  dy  y. 


T  dy 


+  d92  +  dT2  _  -  +  v 
y0  dy  +  yt  dyx  ~  yyx 


x  dy 


(4) 


The  required  trajectories  are  given  in  terms  of  (y^Tj)  or 
in  Appendix  B,  as  are  the  adjoints  required  in  (3).  The  numerical  inte¬ 
gration  of  the  coupled  equations  proceeds  routinely  in  Region  III  of 
parameter  space,  as  there  are  no  unusual  numerical  problems.  It  is  true 
that  de2/dyx  is  nearly  zero  as  02  approaches  6eq  .  but  the  deter¬ 
minant  of  the  coefficient  matrix  of  (4)  is 


+  +  +  + 
V,  ‘  ye\ 


t2[(7v/1-o292/72  +  S02)2  -  Ps02]  >  0  , 
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which  means  that  there  are  no  indetermlnacies  in  the  integration  of  (3) 
and  (4), 


The  gradient  vector  W  in  Region  III  varies  smoothly  with  position 
except  across  the  EDL  ,  and  all  of  the  contours  of  constant  V  inter¬ 
sect  the  capture  circle  outside  the  usable  part.  Note  in  particular  that 
as  the  new  barrier  is  approached,  the  time-to-go  approaches  infinity  con¬ 
tinuously,  while  in  Region  I,  the  time-to-go  is  discontinuous  at  ®  . 


The  safe-contact  trajectories  which  follow  the  tangential  arrival 
at  $2  (see  Fig.  B-l)  are  described  in  polar  coordinates  as 


r  =  -c0  +  yc(ij(-e)  =  0 
e  =  -l  +  (s e+y  \/i-c2e/y2)/0  , 


(5) 


where  r  =  (3  has  been  used  in  the  angular  equation  and  where  0T  <  0 

<  8  .  The  radial  equation  gives  E's  control  for  this  motion  as 

eq 

(cf.  (9)  of  Sec.  3.1): 

\| i=Q  +  cos  ^(ce/y)  ,  (6) 

and,  as  expected,  E's  tangential  angular  velocity  decreases  from  zero 

as  0  decreases  from  0  .  All  tangential  paths  depart  the  capture 

eq 

circle  at  0  =  @T  ,  which  has  been  calculated  as  the  initial  value  of 
02  in  the  EDL  integration.  This  last  portion  of  the  "+"  paths  (see 
Fig.  4.1)  has  E  running  in  a  straight  path  in  real  space,  after  tan¬ 
gential  departure  and  before  termination  at  0  =  0Q  . 


4.3  Conjugate  Point  Replaces  Dispersal  Point 

As  the  speed  ratio  is  increased  for  parameters  above  C  ,  two  ex¬ 
treme  trajectories  from  A  coalesce  to  a  single  curve,  which  is  point- 
tangential  at  @T  >  0up  and  which  ends  at  the  arrival  angle  0=0. 

For  larger  7  ,  such  that  the  y-axis  is  included  in  the  point-tangential 
trajectory,  the  retrogressive  equations  from  the  y-intercept  are 


X  =  1  -  CTj  +  (yi-7T1)st1 

y  =  stj  +  (y1-7T1)ct1  . 


(1) 
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A  conjugal!'  point,  r  ,  is  Jocuted  on  the  iowest  such  trujectory  by  the 
coiulitiou  ss  0  .  Here  it  is  easily  seen  that  the  conjugate 

point  condition  is  met;  i.e«,  when  \ ^  =  y^/y  i  the  Jucobiun  determinant 
is,  using  (1): 


(i-y)st 


<1-7)ct 


=  0 


(2) 


An  K1)L  emanating  from  C  murks  the  right  end  of  trajectories  to  the 
y-axls,  and  the  time  along  these  trajectories  is  <  y  /y  for  those 
above  the  point-tangential  trajectory.  A  trujeetory-siope  discontinuity 
occurs  across  the  EL)L  ,  this  discontinuity  decreasing  to  zero  at  the 
point  C  .  The  tangency  of  the  point-tangential  trajectory  and  the  EDL 
at  C  implies  that  an  indeterminacy  exists  in  the  expressions  for  the 
derivatives  at  tills  point.  This  indeterminacy  is  evaluated  in  Appen¬ 
dix  U. 

The  parameter-space  locus  for  which  the  dispersal  point  A  becomes 
the  conjugate  point  C  is  found  by  requiring  rCp,/,-^)  =  pt  ,  r  CfcA,  7,  ) 

“  0  ,  where  -  (3/ 7  1  using  (i)  with  yA  =  (3  .  Considering  these 

as  two  equations  in  t.t  and  ,  we  determine  us  the  smallest  root 

to 


(l-eij-Jy  TjXy-ctj)  +  (1-7)st  (si  -yt  )  =  0  .  (3) 

When  this  is  solved  numerically  for  1^(7)  ,  the  corresponding  parameter 
is 

t3(7)  =  7Tt  +  (1-7)stt/(7-ct:t)  ,  (A) 

and  a  parameter  locus  is  generated  lor  all  iT  in  the  interval  0  to 
n  .  This  curve  will  be  denoted  us  in  the  parameter  space  drawing 

to  be  shown  as  Eig.  A. 10  in  Sec.  A, 6. 

The  last  point  on  has  the  coordinates  p  =  2  ,  y  =  2/n  =  .635  , 

and  here  C  is  at  the  right  edge  of  the  capture  circle,  with  the  trajec¬ 
tory  departing  tangentially  upward.  This  means  that  1  -  c(y./7)  = 

2  ^ 

3  /2  ,  and  s(y  /y)  =  0  .  Together  these  conditions  imply  that  y^  =  ny 
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and  that  the  conjugate  point  region  is  bounded  from  above  by  the  line, 


(3=2,  this  locus  to  be  denoted  as  in  Fig.  -1 . 10. 

For  parameters  in  Region  111,  the  effect  of  crossing  the  locus  C 
is  shown  in  Fig.  4.2.  The  F1)L  calculation  is  discussed  in  more  detail 
in  Appendix  B,  but  here  it  is  noted  that  the  initial  values  of  y  and 

T2  =  yi/7  '  tT  derlve  ,  rom  «*PandlnB  «•(>!»  XT)  =  Is  ,  ,'(-vi>1t)  =  0  * 
Eliminating  y  ^  leads  to  a  cubic  in  ct^.  , 

2yc3-iT  +  (l+4  7+y2-fi2)c2tT  -  2[l  +  y2+7(l-r32)  ]ctT  +  72(l-p2)  +  1=0,  (5) 
and  the  corresponding  ordinate  is 


>1  =  /xT  +  0-7)8  1^ (r-ctT)  . 
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Comparing  Fig.  4.1  with  Fig.  4.2  illustrates  the  difference  between  tra¬ 
jectories  for  parameters  in  Region  III,  on  either  side  of  the  locus 
given  by  (4). 

4-4  Pursuer's  Dispersal  Line  Replaces  Open  Barrier 

We  recall  that  in  Region  II,  P's  strategy  when  r  =  p  was 
-  +1  with  e  <  9up  and  cp  =  -1  otherwise.  The  barrier  intersection 
with  the  capture  circle  thus  marked  a  discontinuity  in  P’s  strategy. 
For  p  +7  >  1  ,  however,  we  have  seen  in  Sec.  4.1  that  safe-contact 
motion  is  possible  for  cp  =  +1  ,  9  >  0  .  When  parameters  pass  from 

Region  III  to  IV,  the  equilibrium  point  exists  as  long  as  p  <  l  +  y  , 
but  the  barrier  from  this  point  will  be  open,  since  the  parameters  will 
then  fall  above  the  locus  given  by  (4)  of  Sec.  4.1.  This  barrier  starts 
on  the  capture  circle  at  0  ,  and  when  cp  =  +1  ,  the  relative  inward 

motion  will  be  very  slow  when  E  is  located  on  the  capture  circle  just 
ahead  of  0^  .  Indeed,  the  time-to-go  will  approach  infinity,  so  that 
P  should  obviously  turn  away  (<p  =  -1)  when  E  is  close  to  this 
equilibrium  point.  Accordingly,  there  must  be  a  dispersal  line  for  P 
starting  on  the  capture  circle  slightly  in  front  of  the  open  barrier, 
which  then  loses  its  significance. 

Hence,  when  the  barrier  from  0^  is  open,  capture  in  a  finite 

time  is  possible  from  all  initial  conditions,  and  we  may  think  of  this 

barrier  as  being  "replaced"  by  the  PDL  .  We  note  that  the  gradient  and 

P's  switch  function  are  discontinuous  across  this  line,  and  when  the 

equilibrium  point  exists  (i.e.,  jl-J2  <  p  <  1  +  7  ),  the  PDL  begins 

at  an  angle  0  very  slightly  less  than  0  as  given  by  (2)  of  Sec. 
p  eq 

4.1,  When  p  >  1  +  7  ,  so  that  0  no  longer  exists,  the  PDL  may 
begin  on  the  capture  circle  at  an  angle  0p  >  n/2  ,  or  on  the  negative 
y-axis,  or  even  on  the  equivocal  line,  as  will  be  shown.  The  variety  of 
possible  configurations  of  PDL  makes  its  numerical  computation  diffi¬ 
cult,  as  is  further  implied  by  the  following  possibilities: 

i)  The  trajectory  for  cp  =  +1  can  depart  tangentially  from  the 
PDL  ,  indicating  that  the  SE  is  present  for  these  para- 
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meters.  That  is,  the  extension  above  Region  II  of  (see 

Fig.  3.9)  is  required. 

ii)  The  SE  ,  in  turn,  may  be  interrupted  by  the  FL  ;  i.e,,  we 
must  locate  the  extension  of  Cg  of  Fig.  3.9. 

iii)  The  dispersal  point  A  is  replaced  by  the  conjugate  point  C 
for  parameters  to  the  right  of  ,  as  discussed  in  Sec.  4.3. 

iv)  The  point  A  may  separate  the  SL  from  the  EL  ,  or  may 

fall  inside  the  capture  circle,  the  EDL  no  longer  existing. 
That  is,  and  C.  of  Fig.  3.9  must  be  extended  upward 

out  of  Region  II. 

These  considerations  have  an  obvious  significance  in  the  computa¬ 
tion  of  the  PDL  ,  but  for  the  present  we  restrict  attention  to  the 
equations  which  describe  it,  and  to  the  determination  of  those  para¬ 
meters  for  which  it  exists. 

For  parameters  slightly  above  such  that  the  PDL  exists,  the 

configuration  of  exceptional  lines  will  resemble  that  shown  in  Fig. 
4.3(a).  This  configuration  is  representative  of  parameters  in  the  neigh¬ 
borhood  of  (3  =  1.4  ,  7  =  .6  ,  such  that  the  PDL  is  not  interrupted 
by  the  SE  as  mentioned  above.  The  turn-away  region  resembles  that 
which  exists  for  parameters  in  Region  II,  and  when  parameters  are  close 
to  the  locus  C1  ,  the  EDL  and  PDL  are  nearly  coincident,  so  that 
they  can  be  thought  of  as  "replacing"  the  barrier  which  exists  in 
Region  II. 

In  (b)  of  the  figure  is  shown  a  configuration  for  parameters  near 

the  values  8  =  1.82  ,  7  =  .8  ,  so  that  6  does  not  exist,  and  such 

eq 

that  the  point  A  has  been  supplanted  by  the  conjugate  point  C  .  For 
these  parameters,  the  two  dispersal  lines  are  more  widely  separated,  and 
the  turn-away  region  is  proportionately  smaller. 

The  equations  which  specify  the  PDL  are  similar  to  those  described 
in  Sec.  3.3.  That  is,  the  position  is  the  same  for  two  different  retro¬ 
grade  paths,  and  the  main  equation  holds  for  either  of  P's  strategies. 
Using  the  notation  of  Fig.  4.3(a),  for  example,  the  PDL  would  be  found 
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by  simultaneous  solution  of  the  three  equations. 


+  dei  +  dTi 
x9  de2  +  XT  ds2 


X^  +  X  — — 
6  T  d er 


+  d0l  +  dTl  -  -  dT2 

y  -  +  y  —  =  y  +  y  - 

yQ  d e2  T  d02  yG  y t  d02 


v'-v+ 

dy  _  x  x 
dx  V“-V+ 

y~  y 


(a)  p  <  1+7  (b)  P  >  1+7 


FIGURE  4.3 


EDL  and  PDL  in  Region  IV 


At  any  point  x  (g^,T^)  =  x  ,  these  equations  can  be  solved 

for  the  indicated  derivatives,  0^  having  been  taken  as  the  independent 
variable.  The  superscripts  in  (1)  are  chosen  to  agree  with  the  sign  of 
the  switch  function.  Hence,  providing  the  configuration  actually  is  as 
shown  in  the  figure,  it  remains  only  to  find  a  set  of  "initial"  condi¬ 
tions  in  order  for  the  PDL  to  be  numerically  calculable. 

The  first  few  attempts  to  integrate  the  PDL  were  carried  out  for 
parameters  only  slightly  above  ,  so  that  the  relevant  initial  condi¬ 
tions  were  those  at  the  near  end  of  the  PDL  ,  slightly  inward  of  0  . 

eq 

Initial  conditions  could  be  found  only  after  numerical  integration  of 
the  EDL  ,  the  EL  and  the  safe-contact  equations  of  retrograde  motion 
back  to  .  The  initial  conditions  required  for  the  PDL  integra¬ 
tion  would  then  be  x  =  t  =  0  and  9  =  6.  =  9  ,  and  hence  only  9 

1  «  1  &  P  P 

needed  to  be  determined  numerically,  as  described  in  Appendix  B.  This 
procedure  was  followed  by  a  subsequent  PDL  integration,  which  con¬ 
tinued  until  the  switch  function  equalled  zero  on  the  "+"  side  of  it. 

It  was  found  that  the  end  point  of  the  PDL  was  "close"  to  the  end  point 
of  the  EDL  ,  as  pictured  in  Ref.  2.  Subsequent  numerical  studies 
showed,  however,  that  these  endpoints  are  actually  coincident,  and  that 
the  two  dispersal  lines  Join  the  EL  at  a  common  point,  labelled  B  as 
in  the  barrier  case  of  Region  II.  This  means  that  for  configurations 
similar  to  that  shown  in  Fig.  4.3(b),  the  PDL  can  be  profitably  inte¬ 
grated  inward  from  the  far  end,  simultaneously  with  the  EL  .  This  has 
important  numerical  advantages,  particularly  when  p  >  1  +  7  ,  such  that 
the  approximate  location  of  the  near  end  of  the  PDL  is  unknown. 

In  determining  the  parameters  for  which  the  PDL  exists,  we  con¬ 
sider  the  changes  in  its  appearance  due  to  changes  in  p  ,  for  a  fixed 
value  of  7  ,  this  fixed  value  lying  between  0  and  1  .  The  largest 
value  p  for  which  it  can  exist  will  be  found,  and  this  locus  (3(7) 
will  be  denoted  C  ,  for  reference  purposes. 

O 

When  y  is  sufficiently  small,  the  point  A  is  not  a  dispersal 
point,  because  the  retrogressive  switch  time  is 

ts  *  2(n-eo)  <  p/7  *  tA  . 
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Also,  when  the  point-contact  trajectory  is  tangent  to  the  capture-circle 
at  the  angle  t)T  =  Jt/2  ,  the  retrogressive  equations  are 


x  =  1  -  ct  +  (f3-yt)s  (So+t)  =  {3 
y  =  st  +  O-7t)c(0o+t)  =  0  . 


(2) 


But,  at  this  point  we  know  E's  optimal  direction  must  be  tangential, 
ijr  =  9q  +  t  =  it  ,  so  (2)  yields 

X  =  p  =  1  -  ctt  ,  (3) 


where  xT  is  the  least  positive  root  to  the  equation, 

stt  +  ctt  +  nT  ~  1  =  0  •  (4) 

While  this  analysis  appears  unmotivated,  it  may  be  observed  that, 
when  the  independent  variable  is  changed  to  t'  =  T  -  rT  ,  the  trajec¬ 
tory  (2)  is  obviously  symmetric  with  respect  to  the  x-axis.  That  is, 
with  time  t'  measured  from  tangency  at  9Q  =  n/2  ,  and  n  >  the 

trajectory  is  given  by 

X  =  1  +  (P-I)ct’  +  n'BT’ 

(5) 

y  =  -(P-I)st'  +  7t’ct'  . 

This  means  that  the  switch  function  is  also  symmetric  and  is  therefore 
equal  to  zero  at  both  ends  of  the  tangential  trajectory  indicated  in 
Fig.  4.4(a),  This  finally  means  that  parameters  are  such  that  the  PDL 
is  of  zero  length,  having  just  vanished  into  the  capture  circle.  Hence, 
our  sought  function  in  parameter  space  is  obtained  by  eliminating  tt 
from  (3)  and  (4).  In  a  relatively  useful  form,  the  result  is 

7  =  P  -v'ae-P2  (6) 

COS  *  (l-p) 

This  function  denotes  the  left  end  of  the  locus  Cg  (to  be  shown 
in  Fig.  4.10),  and  it  is  appropriate  until  7  is  so  large  that  A  ap¬ 
pears  on  the  capture  circle.  The  criterion  for  this  circumstance  is 

ta  =  2tt  =  p/7 
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(a)  y  v.  .362 


(b)  .362  <*><.  .78 


(c)  .78<r<-  10 


FIGURE  4.4.  Trajectories  for  Parameters  Just  Under  C 


which,  with  (3)  and  (-1),  implies  that  tt  =  cos  1  (-3/5)  =2.21.  It  fol 
lows  that  the  last  point  on  the  locus  (6)  is 


7  =  4/5it  S  .362 
0  =  8/5  =  1.6  . 


(7) 


When  (6)  yields  a  value  of  7  >  .362,  the  point  A  is  outside  the  cap¬ 
ture  circle,  and  is  the  start  of  an  EDL,  The  condition  (6)  must  be  re¬ 
placed  by  another,  when  y  >  .362,  such  that  S  =  0  on  the  inner  side 
of  the  EDL  when  r  =  0  .  That  is,  if  S  =  0  while  r  >  P  ,  a  PDL 
exists  between  this  point  and  the  capture  circle. 

The  conditions  described  above  are  not  subject  to  analytic  treat¬ 
ment,  since  the  EDL  equations  must  be  integrated  numerically.  This 
implies  that  0  is  to  be  adjusted,  at  a  fixed  y  ,  until  the  point  B 
at  the  junction  of  EDL  and  PDL  occurs  at  r  =  p  .  The  qualitative 
appearance  of  the  relevant  curves  for  p  just  below  Cg  is  shown  in 
Fig.  4.4(b).  This  procedure  is  sufficient  as  y  increases  until  the 
dispersal  point  A  is  supplanted  by  the  conjugate  point,  at  y  =  .78  , 
when  C8  is  crossed  by  C?  .  For  larger  7  ,  the  calculations  relate 
to  the  point  C  ,  and  the  conditions  for  the  vanishing  of  the  PDL  are 
slightly  changed.  For  7  greater  than  about  .7,  the  left  end  of  the 
PDL  shrinks  toward  its  right  end  at  the  point  B  ,  as  shown  in  Fig. 
4.4(c).  As  the  speed  ratio  approaches  unity,  it  is  found  that  the  con¬ 
figuration  near  B  is  very  sensitive  to  parameter  changes,  and  conse¬ 
quently  the  precise  location  of  Cg  is  difficult  to  assess  in  this 
region.  The  locus  appears  to  pass  through  the  point  P  =  2  ,  7  =  1 

in  parameter  space,  but  this  observation  has  only  a  numerical  basis. 

We  have  discussed  safe-contact  motion  in  Region  II  for  cp  =  -1  t 
when  P  turns  away  from  E  .  Likewise,  tangential  relative  motion 
appears  for  cp  =  +1  ,  when  E  is  inboard  of  the  equilibrium  point,  for 
parameters  in  Region  III.  For  parameters  below  P  =  1  +  7  in  Region 
IV,  moreover,  safe-contact  motion  can  occur  for  both  of  P's  turn 
directions.  As  P  approaches  Cg  ,  however,  particularly  for  large 
7  ,  the  turn-away  region  does  not  necessary  border  the  capture  circle, 
so  that  tangential  motion  can  occur  only  for  cp  =  +1  .  In  this  case, 
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tangential  motion  appears  whenever  a  point- tangential  trajectory  joins 
the  point  A  ,  or  the  point  C  ,  to  the  capture  circle.  These  trajec¬ 
tories  can  be  expressed  in  terms  of  the  arrival  angle  3o  ,  and  the  retro¬ 
gressive  time-to-tangency  ,  by  expanding  r  =  6  ,  r  =  0  ,  as  two 

equations  in  0*7)  • 

If  the  parameters  are  in  the  region  above  ,  P  never  turns  away 

from  E  ,  and  tangential  motion  can  occur  only  for  cp  =  +1  .  The  maximum 

angle  on  the  capture  circle  for  which  tangential  motion  is  possible  is 

found  by  equating  the  radial  acceleration  to  zero  in  unconstrained  motion, 

with  r  =  3  and  r  =  0  .  Rearranging  the  result  yields  an  equation 

satisfied  by  this  8  >  90°  , 

in 

ri*  =  3s 8  -  (s3  +7 v l-c2e  Z72)2  =  0  .  (8) 

mm  m 

For  parameters  below  C  ,  the  PDL  exists ,  and  the  maximum  angle 

o 

at  which  tangential  motion  occurs  for  cp  =  -1  depends  on  the  subsequent 
trajectory.  We  have  determined  one  such  relation  in  (5)  of  Sec.  3.1  for 
the  case  when  the  departing  trajectory  encounters  the  negative  y-axis. 

In  Region  IV,  the  "last"  tangential  path  for  <p  =  -1  may  encounter  the 
EL  or  it  may  not  exist  at  all.  That  is,  all  paths  departing  the  circle 
with  cp  =  “1  may  do  so  non-tangentially ,  as  shown  in  (a)  of  Fig.  4.5. 

When  the  speed  ratio  is  high  and  (3  >  1  +  7  >  the  shape  of  the  PUL 
is  extremely  sensitive  to  parameter  changes.  Typical  trajectory  varia¬ 
tions  for  7  =  .9  are  shown  in  Fig.  4.5  for  three  slightly  different 
values  of  3  . 


I 

I 


(b)  (3=1.916 


(c)  p=l .918 


FIGURE  4.5.  PDL  Variations  witn  3  for  7  =  .9 


67 


4.5 


Other  Loci  and  Trajectories  in  Regions  IV  and  V 


The  condition  for  the  disappearance  of  the  EDL  (the  extension 
above  Region  II  of  locus  C  )  is  that  the  switch  function  vanish  at  the 
point  A  on  the  point- tangential  trajectory.  Using  the  solution  for 
this  trajectory,  the  condition  requires  satisfaction  of  the  equality 

T.  =  3/7  =  t  =  2  (n-0  )  ,  (1) 

n  SO 

where  0^  and  rT  are  coupled  through  r  =  3  *  f  =  0  ,  as  given  in  (1) 
of  Appendix  B-l. 

The  result  of  eliminating  dQ  and  to  give  3(7)  is  the  seg¬ 
ment  of  C  lying  between  C  and  C  ,  as  will  be  shown  in  Fig.  4.10. 

J  1  O 

In  Fig.  4.6  is  shown  the  nature  of  the  switch  lines  and  trajectories  for 


FIGURE  4.6.  Trajectories  for  Parameters  Between 
C3  and  C5  in  Region  IV 
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parameters  slightly  to  the  left  of  this  portion  of  C  .  Extending  C. 

«3  O 

above  is  accomplished  simply  by  requiring  the  dispersal  point  A  to 

fall  on  the  capture  circle,  using  (1)  of  Sec.  3.2  for  3  <  2  .  For 
(3,7)  left  of  ,  but  below  Cg  ,  the  EL  is  absent  and  the  configura¬ 

tion  is  as  shown  in  Fig.  4.7. 


FIGURE  4.7.  Trajectories  for  Parameters  Left 
of  C5  in  Region  IV 

For  parameters  in  Region  V  above  Cg  ,  the  disappearance  of  the  EDL 
again  results  when  A  is  on  the  capture  circle.  This  is  the  same  func¬ 
tion  as  mentioned  above,  and  it  extends  from  3=7=  0  to  the  point 
P  =  2  ,  7  =  2/n  S  .635  .  It  will  be  seen  in  Fig.  4.10  that  Cg  and  C5 
intersect  C  at  the  same  point  (where  e  =  */2  ). 

o  * 
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The  locus  Cj  denotes  the  upper  border  of  an  area  In  which  the  SE 

exists,  and  it  is  found  by  fixing  y  and  determining  the  largest  p  for 

which  the  required  cotangency  condition  holds.  A  test  of  a  given  pair  of 

parameters  thus  requires  the  retrogressive  integration  of  the  EDL  ,  EL 

and  safe-contact  equations,  until  g  is  found  at  the  base  of  the  PDL 

P 

on  the  capture  circle.  (See  Appendix  B.)  Integration  of  the  PDL  then 
determines  if  the  SE  is  present  by  the  condition  d0  /d^  =  0  .  If 
this  equation  is  not  satisfied  at  any  point  on  the  PDL  ,  the  correspond¬ 
ing  parameters  are  above  the  locus  C 
-  4 

Similarly,  the  locus  Cg  bounds  an  area  of  Region  IV  for  which  the 
SE  is  interrupted  by  a  FL  ,  and  this  locus  is  found  as  in  Sec.  3.4,  by 
Integrating  the  SE  until  the  FL  condition  is  satisfied  or  not.  This 
again  is  a  straightforward  but  tedious  numerical  trial  and  error  pro¬ 
cedure. 

In  Fig.  4.8  is  shown  a  qualitative  sketch  of  the  switch  line  and 
trajectories  resulting  for  a  pair  of  parameters  Just  above  locus  in' 

Region  IV.  Here  it  is  seen  that  the  EDL  meetB  the  y-axis  while  the 
PDL  ,  SE  and  FL  do  not.  This  combination  of  parameters  generates  the 
most  complex  structure  believed  to  exist  in  the  game.  As  many  as  15 
stages  are  possible  for  certain  Initial  conditions,  and  for  these  para¬ 
meters  a  chase  can  include  every  qualitative  feature  of  the  homicidal 
chauffeur  game  except  the  barrier  and  the  switch  line. 

A  brief  discussion  of  the  15-stage  chase  may  be  instructive.  From 
an  initial  position  on  the  PDL  ,  as  shown  in  Fig,  4.8,  P  chooses 
cp  =  +1  until  Just  after  the  forward  portion  of  the  SE  is  encountered. 
The  turn-away  (cp  =  -1)  leads  E  down  to  the  FL  ,  on  which  any  number 
of  switches  by  P  can  be  performed.  This  leads  x  again  onto  the  SE  , 
for  which  cp  =  +1  ,  until  the  PDL  is  met  again,  where  P  now  chooses 
to  turn  away.  Safe-contact  with  cp  =  -1  is  followed  by  a  reverse  chase, 
which  is  continued  along  the  EL  until  the  point  B  is  reached.  Here 
E  chooses  the  UL  strategy  of  the  y-axis  and  therefore  cp  =  +1  until 
this  singular  arc  is  reached,  when  cp  =  0  .  This  leads  x  down  to  the 
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EDL  again,  where  £  chooses  to  run  toward  the  right  equilibrium  point, 
thus  forcing  P  to  take  cp  =  +1  .  The  final  stages  include  a  second 
brief  safe-contact  episode. 


Continuing  our  discussion  of  the  parameter  space  loci  in  Regions  IV 
and  V,  the  locus  continues  to  p  =  2  ,  and  denotes  the  border  of  a 

zone  in  which  a  conjugate  point  exists.  The  array  of  dispersal  and 
equivocal  lines  when  the  point  C  exists,  for  (P,7)  below  Cg  and  to 
the  right  of  C? ,  has  been  shown  in  Fig.  4.3(b). 

The  left-hand  portion  of  Cg  (7  <  .362) ,  above  which  the  PDL  does 
not  exist,  is  given  by  (6)  of  Sec.  4.4.  This  can  also  be  considered  as 
the  locus  of  tangential  contact  at  6T  =  90°  ,  and  it  can  be  extended 
upward  to  p  =  2  ,  as  in  Fig.  B-2 ,  where  it  provides  a  convenient  refer¬ 
ence  line  for  visualizing  the  point-tangent  trajectories. 


We  can  also  find  a  parameter  locus  above  which  optimal  tangential 
motion  does  not  occur  for  «p  -  +1  .  When  7  <  .503  ,  this  locus  is  given 
by  simultaneous  solution  of  three  equations  in  9Q  ,  tT  and  P  (7  be¬ 
ing  a  fixed  parameter).  The  equations  are  r  =  P  ,  f  =  0  ,  and  (8)  of 
Sec.  4.4,  which  is  equivalent  to  ¥  =  0  in  two-dimensional  motion.  This 
locus  will  be  labelled  C1Q  in  the  parameter  space  drawing  of  Fig.  4.10. 
For  7  >  .503  ,  this  curve  is  extended  by  Cg  ,  for  which  A  is  on  the 
capture  circle,  with  tt  =  ta  •  Typical  trajectory  configurations  in 
Region  V  are  shown  in  Fig.  4.9.  Here  it  is  seen  that  P  never  turns 
away  from  E  ,  but  that  safe-contact  and  the  EDL  are  still  possibilities. 
For  sufficiently  large  p  ,  as  in  (d)  of  this  figure,  only  the  y-axis  re¬ 
mains  as  an  exceptional  line. 


4,6  Parameter  Space,  Regions  III.  IV  and  V 

In  completing  our  study  of  the  parameter  space  for  the  game,  it  has 
been  convenient  to  define  three  new  regions.  These  are  bounded  by  the 
heavy  lines  in  Fig.  4.10.  Region  III  features  a  closed  barrier  and  an 
EDL  ,  IV  is  characterized  by  both  PDL  and  EDL  ,  and  for  parameters  in 
V  the  situation  becomes  relatively  simple  again,  the  evader's  dispersal 
line  and  safe-contact  trajectories  disappearing  one  by  one  as  this 
region  is  traversed. 
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Parameter  curve  is  continued  to  7=1,  defining  the  closed 

barrier  parameters  of  Region  III.  Here  the  new  barrier  contacts  the 
capture  circle  at  an  equilibrium  point,  and  in  the  capture  region  there 
exists  an  EDL  which  starts  at  a  dispersal  point,  or  at  a  conjugate 
point,  this  distinction  being  specified  by  a  parameter  locus  . 


FIGURE  4.10.  Regions  in  Parameter  Space 


74 


In  Region  IV  very  complex  games  are  possible,  as  the  barrier  is  re¬ 
placed  by  the  EDL-PDL  combination.  The  upper  boundary  to  this  region, 

C  ,  specifies  the  disappearance  of  P's  turn-away  strategy,  and  the 
8 

other  curves  in  this  region  are  extensions  of  loci  first  defined  for 
Regions  II  and  III.  In  particular,  the  possibility  exists  for  both  the 
SE  and  FL  in  Region  IV,  and  safe-contact  is  optimal  for  both  extreme 
strategies  of  P  . 

The  loci  of  Region  V  are  continuations  of  those  ill  IV,  except  for 
Cg  ,  above  which  the  conjugate  point  disappears,  and  ,  above  which 

safe-contact  for  cp  =  +1  disappears.  The  high  maneuverability  of  P  , 
when  p  is  above  Cg  or  ,  means  that  an  optimal  game  is  simply 

played  by  both  players,  P's  strategy  being  at  most  a  sharp  turn  fol¬ 
lowed  by  a  straight  chase,  and  E's  path  being  straight  throughout. 

The  parameter  char, t  of  Fig.  4.10  also  includes  an  unlabelled  dashed 
line  which  passes  diagonally  through  the  space.  This  locus  corresponds 
to  a  configuration  in  which  the  EL  (or  EDL  )  passes  behind  P  and  to 
the  point  x  =  0  ,  y  =  -p  .  While  not  an  Important  locus  in  the  sense 
of  denoting  the  presence  or  absence  of  switch  lines,  dispersal  lines, 
etc.,  it  does  assist  in  the  sketching  of  qualitative  features  of  a 
specific  game. 
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Chapter  V 


CONCLUSIONS 


It  has  been  shown  that  the  type  of  solution  to  the  homicidal  chauf¬ 
feur  game  is  determined  by  the  values  of  two  independent  parameters. 

The  speed  ratio  and  P's  maneuverability  ratio  are  the  two  parameters 
of  the  game,  and  the  form  of  the  solution  is  believed  to  be  known  for 
all  combinations  of  these  parameters.  An  optimally  played  game  can  con¬ 
sist  of  from  one  to  fifteen  stages,  depending  upon  the  parameters  and 
the  initial  relative  position. 

The  exceptional  lines  which  arise  in  the  game  border  the  regular 
regions  in  the  relative  space,  and  the  exceptional  lines  themselves  may 
or  may  not  be  trajectories.  In  Table  I  are  given  the  distinguishing 
characteristics  of  the  exceptional  lines  which  can  occur  in  the  solution 
to  the  game. 

The  barrier  is  an  optimal  path  when  it  is  closed,  Buch  that  the 
capture  region  is  finite.  When  the  barriers  are  open,  a  range  of  para¬ 
meters  exists  for  which  a  segment  of  the  barrier  is  an  optimal  path, 
along  which  P  lunges  toward  E  .  The  remaining  portion  is  a  locus  of 
initial  conditions  for  which  P’s  strategy  is  a  turn-away.  The  uni¬ 
versal  line,  or  singular  arc,  corresponds  to  straight-line  motion  of 
both  players,  and  is  found  to  occur  for  all  valueB  of  the  parameters. 

Dispersal  lines  for  P  are  found  to  exist  both  as  a  straight  line 
directly  behind  P  and  as  curved  lines  to  the  right  and  left  of  P  . 
For  initial  conditions  on  these  lines,  P  must  choose  between  a  hard 
left  and  a  hard  right  turn,  and  E's  strategy  depends  on  this  choice. 
Similarly,  a  dispersal  line  for  E  is  found  as  a  locus  of  initial  con¬ 
ditions  for  which  E  chooses  from  two  equivalent  strategies.  For  most 
parameters,  P’s  strategy  does  not  depend  on  E's  choice.  However, 
when  a  segment  of  E's  dispersal  line  falls  directly  ahead  of  the 
pursuer,  P's  strategy  does  depend  on  E's  choice. 

Safe  contact  motion,  for  which  E  follows  a  curved  path  in  real 
space,  is  found  to  occur  for  both  extreme  turn  directions  of  P  . 
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TABLE  I 

EXCEPTIONAL  LINES  AND  THEIR  CHARACTERISTICS 


Definition  and  Characteristics 


Barrier,  across  which  W  is  discontinuous.  It  may  be 
open  (Region  II),  in  which  case  it  is  an  optimal  path  only 
if  S  has  the  same  sign  on  either  side  of  it.  It  may  also 
be  closed  (Regions  I  and  III)  and  it  denotes  a  discontinuity 
in  V  in  Regions  I  and  II. 

Universal  line,  a  portion  of  the  y-axis  which  is  a  path  cor¬ 
responding  to  straight  line  motion  of  both  players,  across 
which  W  is  continuous,  and  along  which  S  =  S  =  0  . 

A  portion  of  the  negative  y-axis  which  is  a  pursuer's  dis¬ 
persal  line,  on  which  P  must  choose  between  hard  left  and 
hard  right  turns.  Both  vty  and  S  are  discontinuous 
across  this  line. 

Pursuer's  dispersal  line  for  x  ^  0  ,  which  is  found  by 
numerical  integration.  Both  W  and  S  are  discontinuous, 
and  S  changes  sign,  across  this  line. 

Evader's  dispersal  line  emanating  from  a  dispersal  point  A 
and  across  which  W  and  S  are  discontinuous,  but  S  re¬ 
tains  its  sign  across  it. 

Evader's  dispersal  line  emanating  from  a  conjugate  point  C  . 

Safe  contact  motion  (r  =  (3  ,  r  =  0)  for  which  P  turns 
toward  E  ,  who  follows  a  curved  path  in  real  space  while 
maintaining  safe  contact. 

Safe  contact  motion  for  which  P  turns  away  from  E  . 

Equivocal  line,  along  which  S  =  0  and  E  has  a  choice  of 
two  strategies,  one  of  which  keeps  x  on  the  EL  .  W  and 
S  are  discontinuous  across  this  line. 

Switch  line,  along  which  S  =  0  and  across  which  P's  stra¬ 
tegy  switches.  E's  strategy  across  the  SL  is  continuous, 
since  W  and  S  are  smooth  across  it. 

Switch  envelope,  across  which  W  and  S  are  discontinuous, 
and  which  is  a  trajectory  or  not,  according  to  P's  choice. 

Focal  line,  across  which  W  and  S  are  discontinuous,  and 
which  is  a  trajectory  requiring  either  p  =  +1  or  p  =  -1  , 
according  to  P's  choice. 


Other  exceptional  lines  which  can  also  be  optimal  trajectories  include 
the  equivocal  line,  the  switch  envelope,  and  the  local  line.  While 
these  lines  are  more  difficult  to  explain  or  justify  on  physical  grounds, 
their  appearance  in  the  solution  adds  interest  to  the  problem.  In  fact, 
a  general  conclusion  of  the  study  is  that  a  more  complete  understanding 
of  the  theory  often  follows  from  attempting  a  solution  to  a  particular 
application.  Thus,  for  example,  it  is  unlikely  that  the  switch  envelope 
phenomenon  could  have  been  predicted  from  a  purely  analytical  attack  on 
the  problem. 

As  to  the  parameters  for  which  each  type  of  exceptional  line  occurs, 
Fig.  5.1  labels  twenty  significant  subregions  of  the  parameter  space, 
and  Table  II  summarizes  by  listing  the  exceptional  lines  which  exist  in 
each  subregion,  using  the  notation  given  in  Table  I. 
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FIGURE  5.1.  Definitions  of  Subregions 
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TABLE  II 

DISTRIBUTION  OF  EXCEPTIONAL  LINES 
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C  and  C8  which  separate  the  five  regions  of  the 


parameter 

space 

are  particularly 

significant  because 

i) 

for 

<P,7) 

above 

C1 

,  the  barrier  through  cup  no  longer 

exists,  and 

the  time- 

to-go  is  everywhere  continuous; 

li) 

for 

(P.7) 

below 

S 

,  the  capture  region  is  finite; 

iii) 

for 

0,7) 

above 

C8 

,  P  never  turns  away  from  E  . 

It  is  to  be  noted  that  minor  changes  in  the  configuration  can  occur 
even  in  a  specific  subregion  of  parameter  space.  Thus,  for  example,  the 
configuration  of  Fig.  4.8  occurs  only  for  parameters  in  a  minute  strip 
at  the  lower  edge  of  subregion  IVe.  Similarly,  near  the  upper  edge  of 


79 


subregion  IVf,  the  PDL  may  or  may  not  contact  the  capture  circle.  The 
loci  separating  these  various  possibilities,  however,  are  not  considered 
important  enough  to  be  calculated,  particularly  since  this  would  entail 
very  lengthy  trial  and  error  methods. 
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APPENDIX  A 


SWITCH  ENVELOPE 


The  SE  appears  for  parameters  near  the  lower  boundary  of  Regions  II 
and  IV.  In  Region  II,  it  is  specified  as  beginning  on  the  barrier  at  the 
lower  bound  of  points  for  which  P  saves  time  by  turning  toward  E.  The 
criterion  for  this  occurrence  is 


where 


v+  r  +  v+  r  =  i, 

x  y 


(D 


Vx  =  s(VTo)/D 
V+y  =  c(VTo)/D 


(2) 


♦o  =  eo  '  cos'1(cS0/y) 
d  =Vi-c2eo/r2  (P  +  s eo  - yji-c2eo/y2)  , 

and  where 

x~  =  y  -  7a(0up+tb) 
r  =  1  -  X  -  7c(eup+tb)  . 


(3) 


Retrogressive  trajectories  needed  in  this  study  are  the  barrier, 
constructed  with  q>  =  +1 , 


X  .  I  -  CTb  +  (P  -  7Tb)B(0up+Tb) 

y"  =  sTb  +  (p  “  7Tb)c(VTb)  ’ 

and  the  turn-away  paths  (<p  =  -1)  which  encounter  the  capture  circle 
t"!n;entially , 


(4) 


X+  =  -1  +  CT  +  Ps  (0  -T  )  -  7T  s(f  -T  ) 
o  o  o  o  o  o 

y+  =  ST  +  Pc(0  -T  )  -  7T  c(\|f  -T  )  . 

O  O  O  O  O  O 


(5) 
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The  SE  exists  it  and  only  if  (1)  holds  at  some  point,  Tb,  on  the 

barrier.  To  l'ind  such  a  point,  a  barrier  time  Tb  >  0  is  chosen,  and 

the  resulting  (x,y)  are  solved  for  the  parameters  (0  >*t  )  using  (5). 

These  are  then  used  in  (2)  and  (3)  so  that  the  left  side  of  (1)  can  be 

evaluated.  The  time  t  is  adjusted  until  this  left  side  equals  1,  or 
b 

until  it  is  found  that  the  left  side  is  always  less  than  1,  for  all 
T  <  2(n  -  9  ).  For  parameters  between  and  C4>  the  results  of  the 

iteration  are  initial  values  of  8q,  and  +  =  8up  +  Tb>  which  are 
then  used  as  initial  conditions  in  the  retrogressive  SE  integration. 

The  local  slope  of  the  SE  is  then  expressed  as  a  bilinear  function 


of  the  unknown  Ufi  .'dt  , 
o  o 


dy 

dx 


y  +  y.,(d0  /dt  ) 
J  1  ’Q  o  o 

x  +  x-(d6  /dt  ) 
•t  S  o  o 


(6) 


and  (1)  is  used  with  the  incoming  velocity  on  the  SE, 
x"  =  y  -  7st 


y"  =  1  -  x  -  jef  , 

to  give  an  equation  for  E's  prior  control,  t  .  Following  some  simpli¬ 
fication,  this  equation  is* 

c(f-io+To)  =  =  1  +  f  {cU0-T0>-cVfJ'/1“c2fV72)  •  (7) 

The  prior  control  then  permits  calculation  of  the  SE  slope  in  terms  of 

(0  ,t  ): 

o  o 

dv  _  1  -  x  -  (8) 

dx  “  y  -  7s t 


♦While  multiple-valued^  the  inverse-cosine  operation  is  trouble-free  if 
it  is  recalled  that  i|r  remains  practically  normal  to  the  barrier  as 
the  SE  is  integrated  away  from  it. 


82 


Finally,  (6)  gives 


ddQ  y  -  x  (dy/dx) 

7T  =  "  yQ  -  xfl  (dy/dx)  ' 


and  numerical  Integration  of  (9)  for  a  small  change  in  the  independent 
variable  (here  taken  as  t  )  is  followed  by  updated  values  for  x  on 
the  SE.  The  derivative  in  (9)  is  never  indeterminate.  This  retrogres¬ 
sive  integration  proceeds  until  either  of  two  geometric  circumstances 

arises; 

i)  the  curvature  of  the  incoming  path  equals  that  of  the 
SE,  implying  that  the  SE  is  here  tangentially  extended 
by  a  PDL,  or 

ii)  the  derivative  (9)  equals  zero,  implying  that  the  depart¬ 
ing  trajectory  is  tangential  to  the  SE,  which  in  this 
case  is  tangentially  joined  to  an  FL. 


The  first  case  is  equivalent  to  a  requirement  on  the  time-derivative 
of  i|/  for  motion  along  the  SE,  in  that  the  integration  proceeds  only 
while  d\|T/da#  <  1.  This  derivative  is  expressed  more  conveniently  as 


di|f  _  di|f  M6n  dy 

dTs  “  dy/d0o  dxg  ’ 


(10) 


where,  using  (7), 


djf 

d0 

o 


di|r 

o 

d0 

o 


d  t 


s(+  -VV 


-  [i  + 


■  (♦  +  x  )  dt) 

o  o  o 


d0  =  y0  +  yT  (dTO/dS0) 
o 

=  1  -  x  -  ycf  . 

dT. 

When  dt”/dxg  =  1,  a  PDL  starts,  as  shown  in  Fig.  3.5,  and  it  is  neces¬ 
sary  to  have  the  ad joints  on  the  prior  side,  which  are: 
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(11) 


V  =  Btjr  /D 
x 


V"  =  cf/D 

y 


T„.  denominator  in  (11)  1.  expressed  *■  »' 

on  the  SE.  using  the  main  aquation  on  th.  prior  aid.  to  ilv 


d"  =  (l-x)c'f  +  ya  f  -  7 


(12) 


uni  nt  the  barrier  this  function  vanishes,  and 
Because  the  SE  is  tangential  at  the  barrier, 

tn.  adjoint,  are  lnllntta.  at  this  lose.t  «°1”t  of  «e  se 

The  integration  ol  the  SE  .hen  parametera  are  In  Beglon  IV  1.  much 
th.  ..me,  except  that  the  barrier  no  longer  exlale.  The  1— r  end  ol 
SE  1.  Instead  coincident  s.d  co.ang.nti.l  Pith  .he  upper  end  .<  ■ 

PDL  and  with  a  trajectory  for  <P  =  +1.  aB  "hown  in  Fig‘  4-8’  The  °°" 

dltlo.  denoting  this  point  on  the  PPL  1.  *,/«,  *  » 

d,/.,  .  o>.  The  prior  trajectory  thus  t.ks.  «.  0*“.  ol  th.  b.rrl 

1  w  .  integration  of  the  SE  then  proceeds,  with 

in  the  analysis  above,  and  the  integrate 

no  indeterminacy,  by  calculating  the  derivative  according  to  (9). 
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APPENDIX  B 


COMPUTATION  OF  DISPERSAL  LINES 

B-l  Evader's  Dispersal  Line 

In  the  simplest  case,  which  serves  to  illustrate  the  steps  involved, 
we  may  consider  the  configuration  of  Fig.  B-l,  for  which  the  parameters 
fall  in  Regions  III  or  IV,  to  the  left  of  locus  C? .  The  equations  which 
determine  the  EDL  are  three,  in  the  four  variables  indicated. 

For  a  given  pair  ((3,7)  it  is  first  necessary  to  calculate  the 
initial  values  9^,  i2  for  the  point-tangential  trajectory  leaving  the 
point  A  farthest  to  the  right  in  the  Figure.  This  trajectory  is  read¬ 
ily  expressed  in  terms  of  9 and  the  retrogressive  time  to  tangency 
Tt  such  that  r(0o,TT>  =  £  and  r(0Q,TT)  =  0,  or 


FIGURE  B-l.  Notation  Near  Dispersal  Point  A 
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(1) 


1  -  ctt  +  (t3-7iT)[s(ao+TT)-seo]  -  P7Tt  +  J72t^  =  o 

STX  +  <P-7tT)[  c(0Q+TT)-7l  -  7[B(eo+Tx)-seo1  =  0  . 

Fixing  and  at  convenient  values  permits  calculation  of 

the  corresponding  p,7  shown  in  Fig.  B-2,  to  the  left  of  the  locus  C 
The  initial  value  of  the  time  variable  shown  in  Fig.  B-l  is  Tg  =  P/7-tt- 
When  the  conjugate  point  C  is  present  (£,7  right  of  C?) ,  the  equa¬ 
tions  take  the  form 

1  -  ctt  +  (y1-7tT)stT  +  i[(y1_7V2“  ^  =  ° 

(2) 

(1-7)stt  +  (y1-7'fT)(«.tT  -  7)  =  0  , 

which  can  then  be  solved  for  the  parameters,  given  convenient  levels  of 
yx  and  tT-  These  functions  of  0 , 7)  are  also  shown  in  Fig.  B-2,  to 
the  right  of  locus  C? .  The  points  of  special  interest  in  this  Figure 
are: 

i)  The  tangency  angle  is  0T  =  90°  along  the  diagonal  dashed 
locus,  and  in  general  can  be  found  by 

Pc0T  =  stt  +  (P-7tT)c(0o+tT), 

where  9  and  t_  are  read  from  the  Figure, 
o 

ii)  The  arrival  angle  9  approaches  90°  as  7  approaches  0. 

iii)  The  locus  C1Q  denotes  the  parameters  above  which  tangential 
motion  does  not  occur.  For  7  ^  .503,  determined 

by  the  requirement  that  in  free  motion,  r  =  0  at  tan¬ 
gency,  and  for  .503  <  7  <  .632,  it  is  determined  by  having 
A  fall  on  the  capture  circle.  The  locus  is  completed  with 
the  straight  line  P  =  2,  along  which  tt  =  n. 

In  developing  the  EDL  by  numerical  integration,  it  is  necessary  only 

to  express  the  positions  and  adjointB  as  functions  of  the  parameters 

(y  ,t  )  or  (9  ,t„).  Referring  to  the  general  solutions  of  Chap.  I 
11  2  2 

and  to  Fig.  B-l,  these  are: 
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x+  =  1  -  ct2  +  ps(02+T2)  -  7’12*('l,2+'r2> 
y+  =  st2  +  pc(02+t2)  -  7t2c('l'2+T2) 
x"  =  1  -  ctj  +  (y1-7T1)8't1 

y“  =  stx  +  (y1-7T1)ct1 

Vx  =  s(V'T2*/D 
V*  =  c(^2+t2)/D 

D  ='/l-c2S2/72  (P  -  s92  -  7*/l-c2S2/72) 


(3) 


(4) 


(5) 


V-  =  8t,/(l-7) 

X  1 

v“  =  cij/a-7)  . 


(6) 


These  functions  and  their  derivatives  are  used  in  (3)  and  (4)  of  Sec. 

4.2,  which  are  then  solved  for  the  derivatives  dTj/dy^  dflgA^  and 

dT2/dyl- 

For  parameters  in  the  range  Jl-y  <  P  <  1+7  .  the  equilibrium  point 
0  exists,  and  as  the  EDI.  is  integrated  outward  toward  the  point  B 
(where  it  meets  the  PDL) ,  the  angle  &2  is  very  nearly  equal  to  0eq. 

For  numerical  reasons,  therefore,  a  more  suitable  variable  is  de¬ 
fined  by  e"^  =  0eq  -  02.  as  will  be  discussed  in  App.  B-2. 

In  terms  of  this  variable,  the  slope  of  the  EDL  is  then  express¬ 
ible  as 


%L  =  - 
dx 


D  ST1  -  (1-7)6  ('I'g+tj)  s('leq+T2) 

D  ct1  -  (l-7)c(l2+t2)  C<*eq+V 


since 
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D 


e"^  &  Cdeq/y 


eq 


and  further  aubatltution  permits  cancellation  of  e  from  numerator  and 
denominator  of  the  derivatives;  the  results  are: 


^eq+VV 

dyt  “  "  ?[1  -  C^eq+VT1)] 
df  (yl  ~  nl)c6eq 

<iyl  =  7( i-y)T2(P  -  seeq) 


(7) 


dyl  72t2[l  -  c^eq+VV^ 

The  initial  values  of  the  four  parameters  are  easily  computed  <yj  =>■*  i 
ji=P/7  ,  x2  and  e2  from  Fig.  B-2),  and  the  equations  (7)  are  integrated 

Simultaneously  until  the  switch  function  vanishes  on  the  inner  side  (this 
occurs  long  before  S  =  0  on  the  outer  side)  of  the  EDL,  or  when 

s  =  v>  -  VyX  =  c*eq  -  P»<*oq-«eq>  -  =  0  • 

For  parameters  right  of  C? ,  the  conjugate  point  C  marks  the  near 
end  of  the  EDL,  and  because  the  trajectories  on  either  side  reduce  to  the 
same  path  at  C,  the  necessary  derivatives  are  indeterminate.  Referring 
to  Fig.  B-3,  the  two  paths  give,  at  the  point  C, 

xc  =  1  -  cTl  =  1  -  ct2  +  P s(02+t2 )  -  7t2s(i|r2+T2) 

yc  =  »T1  =  st2  +  Pc(92+t2)  -  7T2c(t2+T2)  , 


and  equating  ad joints  provides  two  more  equations, 
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1-7  =  J l-c2e2/72(p  -  s&2  -  yj  1-c202/72) 

The  last  equation  is  solved  for  >  in  fern's  of  the  parameters, 

and  then  the  preceding  equations  give  initial  values  for 

(1-7)c02 

Tl  =  -r2  +  02  +  cos_1(c  d2/y)  (8) 

yl  =  • 


FIGURE  B-3.  Notation  Near  Conjugate  Point  C 
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The  indeterminacy  at  C  requires  that  second-order  expressions  for 
the  derivatives  be  found.  The  slope  of  the  EDL  is 

djr  t  Vx  -  V* 

^=COt 

y  y 

and  when  expanded  to  second  order  in  At^  At2,  A©2>  we  find 

2  2  2DS 

At?  -  At,  =  — —  A0O  , 

1  2  1  -  7  2  ’ 

where  D0  =  as  given  by  (5).  Likewise,  expressing  changes  in  x 

and  y  to  second-order  yields,  after  some  effort, 

AylA-rl  +  i(1~27)(AT^  -  At2)  =  Oc02/7  -  7t2*0)A02 


AyL  +  (1-7)(Atx  -  Atg)  =  0, 


Eliminating  A02  and  At2  gives  a  large  quadratic  in  ATj/Ay^  but 
when  the  identity 


fJce2/7  -  7T2V9  =  De 


is  used,  the  quadratic  becomes  the  very  simple  linear  expression, 
At 

1  -  2-3 y 


Ayx  2(l-7)(27-l)  * 
The  other  two  derivatives  then  follow  as 


(9) 


At. 


Ayx  “  2(1-7) (27-1) 


(10) 


£0 


Ay 


1  2(27-l)c02[se2  +  y*I  l-c2e0/y‘ 


(1-7)s02 

7(1-c202/72) 


(11) 
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These  expressions  are  useful  only  for  the  first  step  of  the  integra¬ 
tion,  after  which  it  can  proceed  normally,  the  indeterminacy  no  longer 
existing, 

B-2  Pursuer's  Dispersal  Line 

We  first  discuss  the  near  indeterminacy  in  the  construction  of  the 
PDL  which  exists  just  under  the  barrier  for  parameters  near  Cg  in 
Region  II.  For  these  parameters,  it  may  happen  that  (9)  of  Sec.  3.3  is 
"numerically"  indeterminate,  as  it  is  based  on  differencing  nearly  equal 
numbers.  A  straightforward  means  of  overcoming  this  problem  consists  of 
finding  a  point  on  the  PDL  which  is  "far  away"  from  the  indeterminacy. 
That  is,  an  intelligent  guess  is  made  as  to  the  parameters 
which  correspond  to  a  convenient  value  of  0^.  This  value  0Q  is  known 
from  the  Integration  of  the  SE  and  is  taken  back  from  the  end  of  the  SE, 

where  t2  =  0  and  01  =  So '  The  PDL  cont,itions  are  written  at  thi* 
estimated  point  as 


X  *9o’T2*  =  *  (91 'Tl* 

y+<W  =  y(<W  (1) 

v+(WV  =  v'(V 

where  initially  only  0q  is  known.  The  initial  estimates  of  the  other 
variables  allow  calculation  of  the  six  indicated  quantities,  though,  and 
linearization  gives 

x"  60  +  x"  6t.  -  x+  8t  a  x+  -  x“  =  Ax 

0  1  t  1  T  Z 

y“  501  +  y~  6t1  -  y\  St2  a  y+  -  y"  =  Ay  (2) 

60,  +  v"  6t.  -  v+  8t,  a  V+  -  V-  =  AV, 

0  1  T  1  T  Z 

where  the  3x3  matrix  of  partial  derivatives  is  evaluated  using  the 
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current  estimates  of  the  sought  quantities.  Inverting  the  matrix  gives 
corrections  to  the  previous  estimates  as 


-1 


60, 

"  "  + 

X-  X  -X 

—x 

•  w 

Ax 

1 

0  -r  t 

&T1 

= 

y0  yf 

Ay 

.8T2. 

1 

1 

fH 

V 

1 

_ I 

AV 

and  the  steps  are  then  repeated  until  the  calculated  errors  are  negligi¬ 
bly  small.  The  result  of  the  process  Just  described  is  a  set  of  variables 
(0l’Tl’T2^  and  8  corresP°n<lin8  point  (x,y)  on  the  PDL.  The  numerical 
integration  of  the  PDL  can  then  proceed  in  either  direction;  i.e.,  toward 
or  away  from  the  numerical  indeterminacy.  In  a  typical  case,  integration 
toward  the  indeterminacy  is  carried  out  only  until  the  computed  values  of, 
say,  t2  become  noisy  as  the  point  D  of  the  SE  is  approached  from  above. 

We  next  consider  the  numerical  problems  associated  with  the  genera¬ 
tion  of  the  PDL  which  exists  for  (&,y)  in  Region  IV.  For  parameters  in 
the  lower  portion  of  this  region,  the  equilibrium  point  0  is  defined 
on  the  capture  circle,  and  the  near  end  of  the  PDL  begins  just  above  the 
possible  equilibrium  as  in  Fig,  4.3(a).  When  the  parameters  are  closer 
to  the  upper  edge  of  Region  IV,  0  is  no  longer  defined  (according  to 
(2)  of  Sec.  4.1)  and  the  PDL  may  resemble  the  locus  of  Fig.  4.3(b). 


For  the  type  of  configuration  shown  in  (a)  of  this  Figure,  it  is 
qualitatively  clear  that  the  safe-contact  motion  is  very  slow  at  posi¬ 
tions  between  9  and  0_,  because  the  time  on  this  brief  segment  is 
P  T  ^  ^ 

practically  equal  to  the  time  required  for  the  "long  way  around",  via 

the  EL  and  the  positive  y-axis,  say.  Hence,  0  may  be  nearly  equal  to 

0  (where  0  =  0  for  cp  =  +1)  in  which  case  an  appropriate  dependent 
eq 

variable  is  £,  as  defined  in  Section  B-l.  In  the  region  of  parameter 
space  between  (3  =  v  1-y^  and  0  =  1+7,  the  smallest  value  of  this  vari¬ 
able,  £,  is  found  to  be  in  the  range  2  to  5,  corresponding  to  angles 
between  point-tangency  and  equilibrium  of  from  .4°  to  8°,  virtually  at 
the  near  end  of  the  EDL.  At  the  other  extreme,  the  largest  value  of  ' 
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occurs  in  measuring  the  very  small  angle  between  the  PDL  and  the  equili¬ 
brium  point.  A  typical  value  encountered,  for  parameters  near  locus  C  , 

-435  L 

is  £  =  1000,  which  corresponds  to  an  angular  separation  of  10  ra¬ 
dians*.  These  facts  imply  that  £  is  a  necessary  type  of  variable  for 
use  in  numerically  integrating  the  PDL,  when  P  <  1+7. 

The  solutions  to  be  used  for  at  least  the  beginning  of  the  PDL  are 


x+  =  1  —  ctj  +  Ps(01+t1)  -  n1s(+1+‘t1) 
y+  =  StL  +  pc  W1+T1)  - 


(4) 


and 


x“  =  "I  +  ct2  +  MVV  "  yV(W 
y"  =  st2  +  Pc(92“T2)  "  ^2C<V2-T2) 


(5) 


where 

^1  =  91  +  COS 

*2  =  9 2  ”  co»-1<c02/7> 

The  set  (5)  may  require  replacement  by  equations  describing  motion  to  the 
negative  y-axis,  or  to  the  EL,  as  suggested  by  Fig.  4.3. 

According  to  the  definition  of  the  PDL,  we  first  locate  the  angular 

coordinate  9  at  its  near  end,  such  that  the  total  time  to  termination 
P 

is  the  same  by  two  paths.  For  parameters  below  P  =  1+7,  the  configura¬ 
tion  is  as  shown  in  Fig.  B-4. 


*  The  angle  subtended  by  an  electron  at  a  radial  distance  of  one  light- 
-30 

year  is  about  10  radians, 
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Approximating  the  upper  limit  of  integration  in  (6)  by  0  then 

gives 


£  =  -  *n(VV  S  '  ,n<VV  -  W(VVV  - 


(7) 


where 


®eq 


D.(6  )(6-0_ 

6  eq  eq 


)  )  dS 


(8) 


has  been  determined  numerically.  This  procedure  results  in  initial  val¬ 
ues  of  £,  =  T2  =  ®>  and  =  ®p’  88  recluired  for  the  inteRratlon  of 

the  coupled  PDL  equations. 

Because  the  initial  values  of  t[  are  quite  large  (in  contrast  to 
their  values  at  9„) ,  it  is  essential  to  express  both  numerator  and  de- 
nominator  of  the  derivatives  as  factors  of  e  s.  Following  cancellation 
of  this  term,  we  find 

_ V 

d£_  yJ^VT2  a<»ed+V  +  »ty*2Hreq-V  ^  (9) 

dTl  py  ■Jl-c20eq/72  •Jl-c2e2/72  (p  +  S02  -  y  s/l-c202/r2)s(yeq+t1) 


where 


-  ea<6eq+V  -  T'lXV'V  "  V<*eq+Tl>]  ‘ 


In  terms  of  this  derivative,  the  remaining  two  are  written 


Rv  »T  ~  Rx  yT 

\  ye  -  *e 


!Zi_  Rx  ye  ~  Ry  x8 

dTl  x“  yg  -  y™  X0 


(10) 


(11) 
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rw.gital  integration  of  <9)-(ll)  then  generates  the  PDL  until  one  of 
two  things  occurs: 

i)  The  derivative  di^/d^  vanishes,  implying  that  dB^/dtj  =  0, 
so  the  PDL  condition  becomes 


dj,  A  +  ya(d6i/dTi)  yt 

dx  =  \  +  4<dfl1/dti)  "  \ 


Together  with  the  main  equation,  W+l  x+  =  -1 ,  this  implies 
the  SE  condition,  w“-  x+  =  -1.  Hence,  when  d^/d-^  =  0,  a 
switch  envelope  begins.  This  occurs  for  parameters  below  C4 
in  Region  IV. 

il)  The  PDL  meets  the  EDL  at  the  point  B,  where  the  switch  func¬ 
tion  S  =  0.  In  this  case,  for  parameters  above  Cj,  it  is 
often  more  convenient  to  perform  the  integration  of  the  PDL 
from  the  point  B  inward,  rather  than  from  the  capture 
circle  outward. 


In  the  first  case,  the  SE  is  generated  as  described  in  App.  A,  while 
in  the  second  case,  a  different  set  of  equations  is  required  to  generate 
the  PDL.  Referring  to  Fig.  B-5,  we  have  a  relation  between  &x  and  ^ 
from  the  switch  function  equation,  0  =  y  -  V+  x  =  0,  which  yields 


ct^+Tj)  =  7  +  '/l-c201/72  <P  - 


7^1-c^/r2)  .  (12) 


The  object  of  the  following  analysis  is  a  set  of  differential  equations 
relating  *,*„,*.,  and  6  ,  thus  permitting  simultaneous  generation 

®  ..  . _ i _ _ 1 _ _ ,  ..a.4  ahl  A  t.Hfi  PDL 

of  both  EL  and  PDL.  Taking 

equations  give 


i  as  the  independent  variable,  the  PDL 
e 


*  *  x* 

xi  JT  *  \ 

X  *  Hi 

y6  di  +  yt  d-r 

e  e 


+  x 


=  y  +  y 


th 

"2  dTe 


(13) 
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FIGURE  B-5.  Detail  Near  Junction  of  PDL  and  EDL 


and 


where 


dy  = 

dx 


y  +  y  (dT_/dt  ) 
Te  ^2  2  6 

\  +  x‘  (dVdV 


(14) 


Vx  =  ■(♦e"T2>/I>1<x*y) 


V  =  c(t  — r0)/D  (x,y) 
y  e  2  1 

V+  =  s(i|r  +t  )/D(0.  ) 
x  1  1  1 

V*  =  c(t1+T1)/D(01)  . 

At  the  point  B,  of  course,  W  =  W+,  because  *r  =  T_  =  0;  hence 

6  & 

^+*1  =  =  tan  <xB/yB)  ,  and 
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y  •Jl-c2dl/y2)  ■ 


D  ■  -  y/r  -  y  =  D  =  s/l-c  z0  /^(P  +  »©1  - 


The  Indeterminacy  of  (14)  can  then  be  evaluated  by  l'Hopital'a  rule.  Af¬ 
ter  croee-multiplying,  we  find 


-  -  dT2  dSl  dT2  dTl  v  , 

(y-t  +yt  dT)[(V+DxVdr  +  Dx(dr"  +  dr  *  +  >*(cp_x/r  ^ 

e  2  e  e  e  e 

dT2  d8l  dT2  dTl  „  9 

+  (xt  +xt  3t“> £ (Dex-Dy*e)dT-  ‘  Dy(dT  +  37°  "  (rtw)(,,'x/r  )] 


(15) 


=  0 


where  the  subscript  0  i»  an  abbreviation  for  the  partial  derivative, 
d/d©  .  The  pair  of  equations  (13)  is  then  aolved  for 


where 


A1  +  VdVdV 


37  =  A3  +  VdVdV  ■ 
6 


(16) 


=  (y+yr)($-l)/A 


Aj  *  -  2(y+yr)/A 

A3  =  [(r2T1t0-^cei)(l+y/7r)  +  pr  ■s/l-c20l/72($-x/r2)]/A 
A4  =  [(y2i1t0-pc01)(l+y/rr)  -  pr  V l-c^/^d+x/r2)]^ 
A  =  t1[(»61  +  7'/l-cZ01/7i)2  -  PeOJ 


and 


$  =  1  -  7 


r  -x+r 


r(r-y) 
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Substitution  into  (15)  yields  a  quadratic  with  intimidating  co¬ 
efficients, 


dT  2  dx 

<»*>•*■** 


+  C  =  0 


(17) 


where 


A  =  y"  [  (Dey+Dxi|r0)A2  +  Dx(A4+l)]  +  x^  [  (DQx-Dyi|r0)A2-Dy  (A4+l)] 
2  2 

B  =  y“  [  (Dgy+Dxi|r0)A2  +  Dx(A4+l)] 
e 

+  y~  [(D0y+Dx)jra)A1  +  DxAg  +  yx(3$-x/r2)] 


+  x"  [  (D0x-Dyi|f0)A2  -  Dy(A4+l)] 


+  x"  [(D0x-DjT|f0)A1  -  DyAg  -  (r+yy)  (cp-x/r2)] 


and 


C  =  y"  [  (D0y+Dx>jf0)Ai  +  DxAg  +  yx($-x/r2)] 

© 

+  x_  [  (D0x-Dy\|i0 )  A^  -  DyAj  -  (r+yy) ($-x/r2)] 
e 

The  positive  root  to  (17)  is  then  given  by 

dT2/dTe  =  B  -  ^02  "  4AC)/2A>  (18) 

which  permits  calculation  of  the  remaining  derivatives  by  (16) ,  and  the 
local  slope  of  the  PDL  by  (14). 
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APPENDIX  C 


TRAJECTORIES  IN  REAL  SPACE 


It  is  interesting  and  instructive  to  determine  the  motion  of  P  and 
E  in  the  fixed  coordinate  system,  using  (1)  and  (2)  of  Sec.  1.2,  The 
Integration  is  simple  to  carry  out,  requiring  only  <p(t)  and  (f(t)  for 
any  initial  condition,  so  that 


X  (t)  =  Xn  +  /  sen(^  d* 

°J°t 

yp(t>  =  yp0  +f0  C0P(I>  di 


(1) 


where 


0  (t )  =  Sn  + 

P 


o 


<U 


and  t 

xe(t)  =  xe  +  yf  s[^(5)+ep(t>]  di 
°  t 

ye<t)  *  ye  +  7f  cO(8>+0p<S>]  dJ 

o  o 


In  these  equations,  of  course,  q>(t)  and  f(t)  have  been  deter¬ 
mined  by  retrogressive  integration  of  the  equations  of  relative  motion, 
and  either  control  or  both  may  be  discontinuous  when  an  exceptional  line 
is  crossed.  Because  the  relative  trajectory  has  been  computed  for  dis¬ 
crete  time  intervals,  it  is  obviously  Important  to  choose  a  "proper" 
step  size,  such  that  the  intermediate  switch  times  can  be  accurately 
estimated,  and  yet  so  that  the  total  computer  time  required  is  not 
excessive . 

In  real  space,  E’a  motion  is  often  a  sequence  of  straight  lines,  as 
might  be  expected  intuitively.  In  portions  of  the  game,  however,  E  fol¬ 
lows  a  curved  path,  as  in  safe-contact  motion,  and  along  the  EL,  the  SE 
and  the  FL. 

Motion  along  the  EL  (Ref,  1)  has  E  following  a  gradually  curved  pur¬ 
suit  path  in  real  space,  while  P  turns  at  an  intermediate  rate.  Both 
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paths  can  be  interrupted  whenever  E  chooses  to  flee  along  the  tangent  to 
P's  minimum-turn  circle,  whereupon  P  turns  at  maximum  rate. 

The  FL  trajectories  also  present  interesting  real  space  characteris¬ 
tics.  If  P  so  chooses,  his  path  can  be  a  series  of  brief  tangential 
circular  arcs,  while  E's  corresponding  motion  is  an  equal  number  of  near¬ 
ly  straight  Jagged  segments  having  a  sawtooth  appearance.  On  the  other 
hand,  the  SE  trajectory  for  P  is  a  sharp  right  turn,  while  E's  strategy 
thereon  is  a  transcendental  function  of  the  relative  position  which 
amounts  to  a  gentle  curve  in  real  space.  When  P  finally  turns  away, 
causing  E  to  leave  the  SE,  E's  strategy  is  again  discontinuous. 

The  equilibrium  point  in  Region  III  corresponds  to  the  following 
real  space  trajectory:  As  E  traverses  the  barrier  towards  the  equilibrium 
point,  he  is  following  a  straight  path  in  real  space.  After  arrival  at 
the  capture  circle,  however,  with  zero  relative  velocity,  E's  path  is  a 
circle  of  radius  y,  P's  path  in  real  space  being  a  concentric  unit  cir¬ 
cle,  described  at  the  same  angular  rate. 

For  initial  conditions  slightly  inward  of  the  barrier,  E's  path  dur¬ 
ing  safe-contact  motion  is  "nearly"  circular.  That  is,  P’s  motion  is 
circular  and  E  is  moving  very  slowly  tangentially  inward  to  the  usable 
part.  The  same  characteristics  hold  for  parameters  in  Region  IV,  if  E’s 

path  brings  him  to  the  capture  circle  at  an  angle  between  9  and  9  . 

T  P 

In  the  course  of  this  research,  brief  animated  films  were  prepared 
of  multi-stage  chases  for  three  sets  of  parameters,  using  the  computing 
and  plotting  facilities  of  the  Stanford  Computation  Center.  In  the  films, 
the  chauffeur  was  represented  by  the  capture  circle,  and  the  pedestrian 
by  a  small  square.  The  pursuer  was  typically  drawn  in  six  successive  po¬ 
sitions  during  the  time  required  to  traverse  one  diameter  of  the  capture 
circle,  and  each  drawing  was  photographed  four  times,  the  standard  16mm 
projector  speed  being  24  frames/sec.  A  representative  100-second  film 
of  a  chase,  shown  in  both  relative  and  real  coordinate  systems,  required 
some  400  simple  drawings  of  P  and  E,  together  with  a  few  coordinate  lines. 
The  Calcomp  plotting  facility  at  the  Stanford  Computation  Center  was  used 
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for  this  task,  and  subsequent  camera  work  was  ably  executed  by  animation 
specialists  at  the  Stanford  Film  Workshop. 
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